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ABSTRACT 

We compute theoretical predictions for surface brightness fluctuations (SBFs) of single-burst stellar 
populations (SSPs) using models optimized for this purpose. We present results over a wide range of 
ages (from 1 to 17 Gyr) and metallicities (from 1/200 to 2.5 times solar) and for a comprehensive set 
of ground-based and space-based optical and infrared bandpasses. Our models agree well with existing 
SBF observations of Milky Way globular clusters and elliptical galaxies. 

Our results provide refined theoretical calibrations and /c-corrections that are needed to use SBFs as 
standard candles. We suggest that SBF distance measurements can be improved by (1) using a filter 
around 1 fim to minimize the influence of stellar population variations, and (2) using the integrated 
V — K galaxy color instead oi V — Ic to calibrate /-band SBF distances. 

We show that available SBF observations set useful constraints on current population synthesis models, 
and we suggest SBF-based tests for future models. The existing SBF data favor particular choices of 
stellar evolutionary tracks and spectral libraries among the several choices allowed by comparisons based 
on only the integrated properties of galaxies. Also, the tightness of the empirical /-band SBF calibration 
as a function oiV — Ic galaxy color is a useful constraint. It suggests that the model uncertainties in the 
lifetimes of the post-main sequence evolutionary phases are probably less than ±50% and that the initial 
mass function in elliptical galaxies is probably not much steeper than that in the solar neighborhood. 

Finally, we analyze the potential of SBFs for probing unresolved stellar populations in elliptical galax- 
ies. Since SBFs depend on the second moment of the stellar luminosity function, they are sensitive to 
the brightest giant stars and provide complementary information to commonly-used integrated light and 
spectra. In particular, we find that optical/near- infrared SBFs are much more sensitive to the metal- 
licity than the age of a stellar population. Therefore, in combination with age-sensitive observables, 
SBF magnitudes and colors are a valuable complement to metal-line indices to break the age/metallicity 
degeneracy in elliptical galaxy studies. Our preliminary results suggest that the most luminous stellar 
populations of bright galaxies in nearby clusters have roughly solar metallicites and about a factor of 
three spread in age. 

Subject headings: cosmology: distance scale — Galaxy: globular clusters — galaxies: elliptical and 
lenticular, stellar content, formation — stars: late-type, AGB and post-AGB 



1. INTRODUCTION 

When observing the inner regions of a nearby elliptical 
galaxy or the bulge of a nearby spiral galaxy, there are two 
noticeable characteristics of the surface brightness struc- 
ture of these spheroidal stellar systems. The first charac- 
teristic is that the galaxy is brightest in the center with 
the surface brightness falling off gradually with increasing 
radial distance. The second characteristic is only appar- 
ent in good seeing conditions: on small scales, the galaxy 
has a clumpy appearance on the spatial scale of the see- 
ing disk. The dumpiness arises from Poisson statistical 
variations in the number of stars within each resolution 
element. This effect can be easily recognized by visual in- 
spection of images of nearby galaxies like M 31 and M 32, 
in which the small-scale dumpiness can be a few percent 
of the mean surface brightness. Historically, this effect was 
called "incipient resolution." In the modern context, it is 



known as surface brightness fluctuations (SBFs). 

Tonry & Schneider (1988) devised a technique to quan- 
tify SBFs for use as an extragalactic distance indicator 
for undisturbed early-type galaxies. (See also the reviews 
by Jacoby et al. 1992 and Blakeslee et al. 1999.) This 
method relies on using the ratio of the second moment to 
the first moment of the stellar luminosity function (LF) of 
the galaxy as a standard candle: 



(1) 



where rii is the number of stars of type i and luminosity 
Li. The quantity L has units of luminosity and is referred 
to as M when represented as an absolute magnitude. The 
apparent SBF magnitude m can be determined observa- 
tionally, and if the distance to the galaxy is known, M 
can also be determined. Of course, as in the case of or- 
dinary photometry, SBF colors are distance-independent 
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(provided that the fc-correction is neghgible). Using SBFs 
as a distance indicator requires that (1) the bright end of 
the stellar LF in elliptical galaxies and spiral bulges is uni- 
versal, or (2) variations in the LF from galaxy to galaxy 
can be measured and corrected so that M remains a stan- 
dard candle. 

SBFs are an intrinsic property of a stellar population 

as a whole. Therefore, in addition to their utility as a 
distance indicator, SBFs offer much promise in adding to 
our knowledge of the stellar content of elliptical galaxies. 
In fact, the use of SBFs for stellar population studies ar- 
guably preceded its use as a distance indicator: Baum 
& Schwarzschild (1955) used the " coimt-brightness ratio," 
the ratio of the number of resolved stars to the integrated 
light, to study the populations of M 31 and M 32. Hence, 
the idea of using observations near or at the limit of reso- 
lution to explore stellar populations has a long history. 

Furthermore, SBFs provide information about stellar 
populations unique from ordinary integrated light, which is 
the first moment of the stellar LF. Since SBFs also depend 
on the second moment of the stellar luminosity function, 
they are especially sensitive to the most luminous stars 
in elliptical galaxies, the evolved cool giant stars. Thus, 
SBFs can put stronger constraints on the evolution of these 
stars than integrated light alone. Both the interior struc- 
ture and emergent spectral energy distributions of cool 
giant stars are poorly understood, especially near the tip 
of the red giant branch (RGB) and the asymptotic giant 
branch (AGB) populated by low and intermediate mass 
(< 5 — 7Mq) stars. Ideally, we would study these stars us- 
ing Local Group star clusters, which comprise populations 
of homogenous composition and age. However, because 
cool giants evolve rapidly, only a handful are present in any 
cluster; therefore, small number statistics and stochastic 
fluctuations are undesireable factors (e.g., Santos & Fro- 
gel 1997). In this context, SBF analyses of entire galaxies 
can complement star cluster studies, since galaxian light 
arises from several orders of magnitude more stars. Also, 
while nearby globular clusters arc mostly metal-poor, the 
dominant stellar populations in ellipticals are thought to 
be generally old and metal-rich. There is a dearth of such 
systems in the Local Group — the best examples are the 
bulges of the Milky Way and M 31, but these may be im- 
perfect analogs. In order to study metal-rich stellar evo- 
lution, one naturally turns to elliptical galaxies. 

There are two basic motivations for modeling SBFs: (1) 
as pointed out by Tonry & Schneider (1988). one can de- 
rive the calibration of SBF absolute magnitudes purely 
from models provided that stellar populations in galaxies 
can be modeled accurately, and alternatively, (2) one can 
use the observations of SBF magnitudes and colors to test 
and improve the models. 

The first attempt at deriving a purely theoretical SBF 
zeropoint was that of Tonry et al. (1990). They used stel- 
lar evolutionary models from the Revised Yale Isochrones 
(Green et al. 1987, hereinafter RYI) supplemented with 
simple prescriptions for the horizontal branch and AGB; 
the resulting zeropoint and especially its dependence 
on the integrated V — Ic galaxy color disagreed signifi- 
cantly with observations (Tonry 1991). This was due to 
the fact that the RYI giant branches failed to turn over 
in the optical at high metallicity, probably because of in- 



accurate bolometric corrections for the coolest giant stars 
(Mould 1992; Ajhar & Tonry 1994). A subsequent study 
was made by Worthey (1993a), who computed SBF magni- 
tudes using his own population synthesis models (Worthey 
1994). For the main-sequence and RGB stars, these mod- 
els used an amalgamation of isochrones from Vanden- 
berg and collaborators with the RYI; post-RGB evolution 
was added using the fuel consumption theorem, including 
"schematic" treatments of the HB as a single red clump 
and of AGB evolution using a variety of theoretical pre- 
scriptions. The resulting SBF predictions agreed well with 
the observed optical SBF colors {V, Re, and Ic) and with 
the empirical calibration of Ic versus V — Ic (see also Tonry 
et al. 1997). Buzzoni (1993) also computed predictions for 
SBF magnitudes which were consistent with the existing 
optical SBF data at the time, though there were some un- 
certainties in transforming from the Johnson filters used 
in his models to the Kron-Cousins ones used for the ob- 
servations. Since his models used older theoretical spectra 
(Bell & Gustafsson 1978), Buzzoni had to extrapolate the 
spectra for wavelengths longward of 1.08 /im and also for 
stars with Tg/y < 4000 K. For both of these reasons, the 
Buzzoni (1993) models are expected to be less accurate for 
SBF predictions in the IR. For example, their predictions 
are at least several tenths of a magnitude fainter in the 
/T-band than the observations. 

The major observational effort on SBFs has been fo- 
cused on /-band measurements for distance determina- 
tions. As mentioned above, comparisons of stellar popula- 
tion models with observed SBF magnitudes and colors can 
also help us calibrate the colors (i.e., stellar spectral en- 
ergy distributions [SEDs]) and numbers (i.e., evolutionary 
lifetimes) of cool luminous giant stars in old stellar popu- 
lations. SBF stellar population studies have been less ex- 
plored than distance measurements, partly because of the 
lack of suitable datasets. Multicolor optical (VRI) SBF 
measurements for Virgo cluster galaxies from Tonry et al. 
(1990) were analyzed by both Worthey (1993a) and Buz- 
zoni (1993), who reached opposite conclusions on whether 
the optical SBF colors indicated that the galaxies con- 
tained a significant metal-poor component. Ajhar & Tonry 
(1994) measured /-band SBFs for a sample of Galactic 
globular clusters, the only observations to date for these 
systems. They sought to understand the empirical /-band 
SBF zeropoint and its correlation with integrated V — Ic 
galaxy colors, as well as the conflict between the data and 
the RYI models; they also addressed the use of optical 
SBF colors for stellar population studies in galaxies and 
for disentangling age from metallicity effects in globular 
clusters. Sodemann & Thomsen (1995, 1996) performed 
a detailed study of the optical SBF gradients within the 
nearby ellipticals M 32 and NGC 3379. Finally, Jensen 
et al. (1998) found reasonable agreement between their 
/f '-band SBF data for 11 nearby galaxies and the Worthey 
(1993a) models, with most of the galaxies lying around the 
[Fe/H] = -0.25 models with a spread in ages. 

Now is a ripe opportunity to revisit the issue of stel- 
lar population modeling of surface brightness fluctuations. 
There have been significant recent improvements in the 
stellar evolution calculations and spectral libraries used by 
population synthesis models. For example, the latest stel- 
lar evolution calculations include updated input radiative 
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opacities (e.g., Iglesias et al. 1992). There has been even 
more progress on the observational front. The amount 
of 7-band SBF data has increased by nearly tenfold since 
the early modeling of Worthey (1993a), and new data in 
the near-infrared have extended the spectral range of SBF 
measurements (Pahre & Mould 1994; Jensen et al. 1998; 
Liu et al. 2000; Mei et al. 2000). 

In this paper, we present new models for opti- 
cal/infrared SBFs of intermediate-age and old single-burst 
stellar populations (SSPs) and discuss their implications 
for SBF distance measurements and stellar population 
studies. Though interesting issues remain to be addressed 
by blue/near- UV SBF measurements (Worthey 1993b). we 
focus on the optical and near-infrared (i.e., V-band to K- 
band) SBFs, since these constitute; the bulk of past and 
ongoing observations. In § 2, we describe our models for 
computing SBFs. The models cover a wider range of ages 
and metallicitics than in previous SBF studies. Further- 
more, we have optimized the models for this work by re- 
fining the prescription for the luminous cool stars, which 
are important contributors to the SBF signal. In § 3, we 
present the predictions of our models, including SBF mag- 
nitudes, integrated colors, and the fractional contribution 
of difi'erent stellar evolutionary phases to the SBFs. We 
also derive theoretical calibrations and fc-corrections that 
are needed to use SBFs standard candles. In § 4, we com- 
pare our results with current SBF observations and discuss 
implications for the stellar content of elliptical galaxies. In 
§ 5, we review the uncertainties in our results and explore 
the potential of SBFs for breaking the age-metallicity de- 
generacy in studies of elliptical galaxies. Finally in § 6, we 
summarize our findings and offer some future directions 
for SBF studies. 

2. SYNTHESIZING SBFS OF STELLAR POPULATIONS 

In this section, we present our models for computing 
SBFs. We first recall the origin of the SBF signal of stel- 
lar populations. We then review current observations of 
the stars which dominate this signal in order to establish 
a framework to interpret the model results. Readers inter- 
ested only in the description of our models should skip to 
§ 2.3. 

2.1. Origin of the SBF Signal 

Since SBFs are weighted by the square of the stellar 
luminosity, they are very sensitive probes of the most lu- 
minous cool giant stars. For example, for the idealized case 
of a stellar population with a simple power-law luminos- 
ity function, one can easily show via equation (1) that the 
fluctuation huninosity scales linearly with the maximum 
luminosity of the stars. Figure 1 illustrates the greater 
sensitivity to luminous giant stars of SBFs compared to 
ordinary integrated light for a more realistic model stel- 
lar population (see § 2.3 below). For this 12 Gyr solar- 
metallicity model, about 90% of the SBF signal is con- 
tributed by the brightest 2 mags of the stellar luminosity 
function in the near-IR and the brightest 3-5 mags in the 
optical. Clearly, the optical/infrared SBF signal is heavily 
weighted to the very brightest stars in a stellar population. 

^Observations actually measure the brightest RGB stars in the clusters, rather than the true tip of the RGB. Therefore, there is a potential 
systematic underestimate of the tip luminosity due to statistical sampling. Unpublished estimates by Rood & Crocker (1997) show that this 
effect is probably not significant for the Progel et al. data (see also Ferraro et al. 1999). 



2.2. Observed Properties of Cool Giant Stars 

Since SBF magnitudes and colors are expected to closely 
follow the peak magnitude^ and colors of the giant branch, 
it is insightful to review the observed properties of cool 
giant stars, both on the RGB (H-shell burning) and the 
AGB (He-shell burning), as a function of age and metal- 
licity for Milky Way and Magellanic Cloud star clusters. 

2.2.1. Giant Branches of Old Stellar Populations 

The colors of the giant branches of old (>5 Gyr) stellar 

populations are primarily driven by metallicity, with more 
metal-rich systems having redder giant branch colors (e.g., 
Frogel et al. 1983). This is illustrated in Figure 2, which 
shows the change in the temperature of the RGB and AGB 
for a change in metallicity from [Fe/H] = -0.7 to -1-0.4 
with ages of 5 and 12 Gyr. This trend results from the in- 
crease in opacity with metallicity from ions, which are 
the principal source of continuum opacity for stars with 
Teff « 3000 — 6000 K. Metals with low-ionization poten- 
tial are the primary electron donors. Thus, as the metal- 
licity increases, H" opacity increases and the giant branch 
temperature drops. It is interesting to note that since Mg 
and Si are also significant donors along with Fe (Renzini 
1977), the giant branch color traces the total metallicity 
[Z/H] and not just the Fe abundance [Fe/H] (e.g., Geisler 
1984; see discussion in Salaris & Cassisi 1996). In stellar 
evolution models, the absolute temperature of the RGB 
depends on the choice of the mixing length used to param- 
eterize the interior convection (e.g., Renzini & Fusi Pecci 
1988). 

Similarly, the shape of the RGB, as measured from the 
RGB slope in optical/IR color-magnitude diagrams, is also 
primarily driven by the metal abundance for old popula- 
tions, with more metal-rich clusters having steeper slopes 
(Kuchinski et al. 1995; Kuchinski & Frogel 1995; Tiede 
et al. 1997; Ferraro et al. 1999). This implies that giant 
stars which contribute to the SBF signal will have a larger 
spread in color for metal-rich populations than metal-poor 
ones. 

The tip of the RGB (TRGB) is delineated by core-He 
ignition. The TRGB bolometric luminosity in Milky Way 
globular clusters increases modestly with metallicity (Fro- 
gel et al. 1981, 1983).^ The optical and IR magnitudes for 
the TRGB have opposite dependences on metallicity. In 
the optical, increasing metallicity leads to increased opac- 
ity from molecular lines, especially from TiO, which re- 
duces the optical flux. Hence, the TRGB in the optical 
becomes cooler and fainter as metallicity increases, first 
turning over in the V-band and then also in the /-band 
for the most metal- rich Milky Way globular clusters (e.g., 
Lloyd Evans & Menzies 1977; Da Costa & Armandroff 
1990; Ortolani et al. 1991; Bica et al. 1991). Thus the 
brightest stars of metal-rich globular clusters are fainter 
in the optical than those of metal-poor clusters. On the 
other hand, in the near-IR the i^-band magnitude of the 
TRGB of Milky Way globular clusters rises monotonically 
with metallicity, with a roughly 1 mag increase over the 
range -2.2 < [Fe/H] < -0.2 (Ferraro et al. 1999). 
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The AGB is expected to follow the behavior of the RGB 
in old populations and reach cooler temperatures with in- 
creasing metal abundance. AGB studies in globular clus- 
ters are hampered by the scarcity of stars found in this 
short evolutionary stage in any given cluster. The AGB is 
usually divided into two phases: (1) the initial early AGB 
(E-AGB), where He-shell burning proceeds steadily, and 
(2) the thermally pulsing AGB (TP- AGB), characterized 
by oscillations in luminosity due to periodic flashes of the 
He shell. Mass loss is believed to occur during both of these 
phases and to end in a superwind stage. In metal-poor 
([Fe/H] < —1) Milky Way globular clusters, the luminosi- 
ties of the brightest AGB stars do not exceed the tip of the 
RGB. In more metal-rich systems, long-period variables 
are observed above the tip of the RGB, which arc presum- 
ably TP- AGB stars (Frogel & Elias 1988). In fact, since 
AGB stars follow a core mass-luminosity relation (Paczyn- 
ski 1971) and the main-sequence turnoff mass increases 
with metallicity at a fixed age (Renzini & Fusi Pecci 1988), 
higher metallicity leads to brighter AGB stars. 

2.2.2. Extended Giant Branches of Intermediate-Age 

Populations 

Information on intermediate-age (« 0.5 — 5 Gyr) stellar 
populations largely comes from studies of star clusters in 
the Large and Small Magellanic Clouds (LMC/SMC). The 
brightest stars in these clusters are carbon AGB stars with 
bolomctric magnitudes much brighter than the tip of the 
RGB (Mould & Aaronson 1979; Frogel et al. 1980; Mould 
& Aaronson 1980; Cohen et al. 1981). These very cool 
stars provide a significant fraction of the total luminos- 
ity of the clusters (Pcrsson et al. 1983). They are redder 
and more luminous than the M-type AGB stars. In con- 
trast, the brightest stars in old LMC/SMC clusters are 
M giants no brighter than the brightest stars in Galactic 
globular clusters (Aaronson & Mould 1985; Frogel et al. 
1990; Marigo et al. 1996). As a result, the near-IR colors 
of intermediate-age clusters (types IV to VI of Searle et al. 
1980) are much redder than those of older or younger clus- 
ters, whereas the optical colors redden monotonically with 
cluster age. 

Carbon stars are believed to arise from carbon dredge- 
up during the TP- AGB phase (Iben & Renzini 1983). The 
observed ratio of C-stars to M-stars in Local Group galax- 
ies anti-correlates with metallicity, C-stars being the dom- 
inant spectral type for [Fc/H] < —1 (Blanco et al. 1978; 
Cook et al. 1986; Pritchet et al. 1987). This occurs because 
stars with lower metallicities require fewer thermal pulses 
to change from an 0-rich to a C-rich star. In addition, 
dredge-up is more efficient at lower metallicities because it 
begins at lower core masses, i.e., at earlier ages. There are 
no examples of solar-metallicity intermediate-age clusters 
in the LMC/SMC which would be most directly compara- 
ble to elliptical galaxies. However, models by Vassiliadis & 
Wood (1993) predict that the bolometric luminosity of the 
tip of the AGB should remain much brighter than the tip 

^Extended giant branches have been claimed to be detected in optical and near-IR images of M 32 and the bulge of M 31 (Rich & Mould 
1991; Preedman 1992; Elston & Silva 1992; Rich et al. 1993), which would mean recent star formation in these supposedly canonical examples 
of old stellar populations. However, more recent work has suggested this conclusion is erroneous and due to severe image crowding (Depoy 
et al. 1993; Grillmair et al. 1996; Renzini 1998; Sodemann & Thomsen 1998; Jablonka et al. 1999). In addition, as pointed out by Renzini 
(1998), the presence of objects brighter than the tip of the RGB {M^gi « —4 for solar-metallicity) cannot be readily interpreted as the sign of 
an intermediate-age population. Long-period variables in metal-rich globular clusters do reach Mbol ^ ~5i so the only unambiguous signs of 
an intermediate-age population would be either the presence of (1) stars with even larger bolometric luminosities or (2) many more stars with 
^bol = ~4 to —5 than expected in a metal-rich population. 



of the RGB and have even cooler temperatures and redder 
colors than in LMC/SMC clusters (see also discussion in 
Silva & Bothun 1998). 

Since the giant branches of intermediate-age Magellanic 
Cloud clusters reach higher luminosities than those of 
old Galactic clusters, they are often referred to as "ex- 
tended" giant branches and are considered the hallmark of 
intermediate-age populations (e.g., Guarnieri et al. 1997). 
This occurs because the core mass-luminosity relation for 
AGB stars implies that younger poptdations have more lu- 
minous AGB stars. In the context of SBF studies, ellipti- 
cal galaxies are generally thought to be composed entirely 
of old stellar populations. However, any recent episodes 
of star formation would lead to extended giant branches. 
The observable consequence of this would be brighter fluc- 
tuation magnitudes (and perhaps somewhat redder fluctu- 
ation colors), especially in the infrared, without significant 
changes in the integrated colors of the galaxies.^ We re- 
turn to this point in § 4.3 and § 6. 

2.3. Stellar Population Synthesis Models for SBFs 

We use the latest version of the Bruzual- Chariot popula- 
tion synthesis models (Bruzual & Chariot 2000, hereinafter 
BC2000). These models allow us to predict the spectral 
evolution of stellar populations with arbitrary star forma- 
tion rates and initial mass functions (IMFs) for a wide 
range of metallicities. The main ingredients of the models 
are the stellar evolution theory used to predict the dis- 
tribution of stars in the theoretical Hertzsprung-Russell 
diagram and the library of spectra assigned to stars as 
a function of effective temperature and luminosity. The 
integrated spectrophotometric properties of an entire stel- 
lar population are obtained by summing the spectra of 
its component stars. We now describe the model features 
most relevant to the present work. 

2.3.1. Model Inputs 

The BC2000 models present several major improve- 
ments over the earlier models of Bruzual & Chariot (1993). 
They also offer several choices of input stellar evolution 
theory and atmospheres (Table 1). Stellar evolution can 
be followed according to the prescription of the Geneva 
school for solar metallicity or that of the Padova school 
for metallicites ranging from 1/200 to 2.5 times solar. 
(See Chariot et al. 1996 for a detailed description of the 
two prescriptions, and Bruzual et al. 1997 for a complete 
list of references.) These prescriptions include all phases 
of stellar evolution from the zero-age main sequence to 
the beginning of the TP- AGB (for low- and intermediate- 
mass stars) and core-carbon ignition (for massive stars). 
In both sets of tracks, solar metallicity corresponds to 
Zq = 0.020, and in the Padova tracks, the metallic- 
ity scales as [Fe/H] = 1.024 \ogZ + 1.739 (BerteUi et al. 
1994). We note that the Padova and Geneva models pre- 
dict very different fractional contributions from RGB and 
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AGB stars to the integrated light of solar-metallicity SSPs 
(Bruzual 1996); this has strong implications for SBF pre- 
dictions.'* 

Bruzual & Chariot (2000) supplemented these tracks 
with a new prescription for the thermally-pulsing regime 
at the tip of the AGB, which improves over that adopted 
by Bruzual & Chariot (1993). This prescription was mo- 
tivated by the inability of that proposed by Chariot & 
Bruzual (1991) to account for the SBF properties of ob- 
served galaxies when we started the prescint study. In 
the new prescription, the effective temperatures, bolomet- 
ric luminosities, and lifetimes (of both the optically-visible 
and superwind phases) of TP- AGB stars of various metal- 
licities are taken from the models of Vassiliadis & Wood 

(1993) . For stars undergoing a transition from oxygen- 
rich (M-type) to carbon-rich (C-type) on the TP- AGB, the 
relative lifetimes of the two phases as a function of metal- 
licity arc taken from the models of Groenewegen & de 
Jong (1993) and Groenewegen et al. (1995). These mod- 
els reproduce the observed ratios of C to M stars in the 
LMC and the Galaxy; however, since they do not extend 
to sub-Magellanic (Z < 0.004) nor super-solar {Z > 0.2) 
metallicities, the same relative durations of the M-type 
and C-type phases are applied to the Padova tracks at the 
more extreme metallicities. Finally, the tracks are supple- 
mented with a prescription for post- AGB evolution and 
with unevolving main-sequence stars in the mass range 
0.1-0.6 Mq. 

The BC2000 models also offer a choice of spectral li- 
braries, summarized in Table 1. For solar metallicity, a 
quasi-empirical set of spectra is available from a compila- 
tion by Pickles (1998). The spectra of MO -MIO giant stars 
are the only non-empirical ones in this library, as they are 
based on the synthetic M-giant spectra computed by Fluks 
et al. (1994). These models were constructed to reproduce 
period-averaged spectra from observations of long-period 
variable TP- AGB stars. In the remainder of this paper, 
we refer to this spectral library as the empirical SEDs. 

For stars of all metallicities, a comprehensive spectral 
library has been assembled by Lejeune et al. (1997, 1998, 
hereinafter LCB97). It is a combination of model atmo- 
spheres by Kurucz (1995; private communication to Leje- 
mic), Bessell et al. (1989, 1991), and the same Fluks et al. 

(1994) spectra as above, rebinned onto homogeneous scales 
in wavelength and physical parameters (effective temper- 
ature, surface gravity, and metallicity). Hereinafter, we 
refer to this spectral library as the theoretical SEDs. 

LCB97 further corrected the synthetic spectra in their 
library by requiring that these reproduce the color- 
temperature relations observed for solar-metallicity stars 
(see also Lejeune c!t al. 1996). This adjustment is espe- 
cially important for the M-star spectra. Lacking empiri- 
cal calibration for non-solar metallicities, LCB97 applied 
the corrections derived at solar metallicities to all other 
metallicities. We refer to this spectral library as the semi- 



empirical SEDs. 

Bruzual & Chariot (2000) supplemented all three of 
these spectral libraries with period-averaged spectra for 
C-type stars on the TP- AGB based on model atmospheres 
from Loidl et al. (1999, private communication; see Hofner 
et al. 2000). BC2000 applied empirical corrections to these 
spectra to match observed color-color calibrations by Men- 
doza & Johnson (1965). The spectra of stars in the super- 
wind phase at the end of the TP- AGB are based on obser- 
vations by Lc Sidaner & Le Bertre (1996) and Le Bertre 
(1997, and references therein). 

The models used here have been tested successfully 
against observed spectra of star clusters and galaxies 
(Bruzual et al. 1997; Bruzual & Chariot 2000). Comple- 
mentary descriptions and applications of previous versions 
of these models can also be found in Chariot et al. (1996), 
Bruzual et al. (1997), and Kauffmann & Chariot (1998). 
Bruzual (1996) has compared the results of different pop- 
ulation synthesis codes and examined the dependence of 
integrated galaxy colors, SEDs, and spectral indices pre- 
dicted by the Bruzual & Chariot (1998) models (very sim- 
ilar to BC2000 for his purposes) on the choice of input 
evolutionary tracks and spectra. 

Unless otherwise indicated, in all calculations below we 
adopt a Salpeter (1955) IMF truncated at 0.1 and 125 Mq. 
For old (>5 Gyr) populations, adopting the Scale (1986) 
IMF would lead to fluctuation magnitudes only <0.05 mag 
brighter and integrated colors only <0.02 mag bluer. (The 
changes are somewhat larger at younger ages; at 1 Gyr, 
the Scalo IMF leads to M's fainter by about 0.1-0.2 mag.) 
These changes would have a negligible effect on our results 
and the agreement of our models with SBF observations. 
(See § 5.2 for a discussion of alternative IMFs.) Further- 
more, in this paper, we approximate star clusters and el- 
liptical galaxies as single-burst stellar populations. While 
this approximation is adequate for our present purposes, it 
limits investigations of the star formation history of ellip- 
tical galaxies. This is our goal in a complementary study 
(M. Liu et al., in preparation; see also § 6). 

2.3.2. Model Outputs 

We compute SBF magnitudes of single-burst stellar 
populations for ages of 1 17 Gyr and metallicities Z = 
0.0001-0.05 {Z/Zq = 0.005-2.5, [Fe/H] -2.4 to +0.4). 
We generate results for the BVRJcJHK'KsKLL'M 
bandpasses used by ground-based observatories and the 
F814W^, FllOM, FllQW, F160W^, and i^222M bandpasses 
used by the WFPC2 and NICMOS instruments aboard 
HST. The zeropoints for our magnitudes are computed on 
the Vega magnitude system. 

The response functions for the Johnson-Morgan BV 
bandpasses are from Buser (1978). The Rc and Ic band- 
passes are from Bessell (1990). For the ground-based 
near-infrared JHKLL'M bandpasses. we use traces from 
Bessell & Brett (1988). These arc similiar to the filter sys- 

^Girardi et al. (2000) have recently updated the Padova library. The new library covers only a limited range of initial stellar masses and 
metallicities and cannot be combined with the library used here because of the slightly different chemical composition used at a fixed metallicity. 
The new tracks are almost identical to the previous ones, the main purpose of the new library being a finer sampling of initial stellar masses 
(Girardi et al. 2000; A. Bressan 2000, private communication). 

^Preliminary results for the iiT-band SBF magnitudes in Liu et al. (1999a,b) used the old TP-AGB prescription and are supplanted by the 
results in this paper. 

®There is potential ambiguity in the notation of L and M. In the case of L, this can be either the fluctuation luminosity or the SBF 
magnitude in the L-band filter {Xeff = 3.5 /im); we discuss the filter only in § 5.3 and Table 2. In the case of M, we use this to refer 
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tern used at the UK Infrared Telescope (UKIRT), which 
is described in Leggett (1992) and Casali & Hawarden 
(1992). We choose not to use the CIT bandpasses (Elias 
et al. 1982) as the effective wavelength of the CIT J-band 
is noticeably redder than that of most other IR systems 
(e.g., Leggett 1992; Bessell & Brett 1988) due to a Si lens 
in the original CIT dewar which defined the blue end of 
the bandpass. 

We also include calculations for popular alternatives to 
the K filter: the K' (1.9-2.3 micron; Wainscoat & Cowic 
1992) and (2.0-2.3 /xm; McLeod et al. 1995) band- 
passes. The Kg filter trace is from the IR imager CIRIM 
used at Cerro Tololo Interamerican Observatory, and the 
K' filter is from the UCLA Gemini instrument at Lick Ob- 
servatory (McLean ct al. 1994). Our K' filter is virtually 
identical to that in the Hawaii QUIRC infrared camera 
(Hodapp et al. 1996) used by Jensen et al. (1998, 1999) 
for their SBF measurements. However, it is slightly differ- 
ent than the trace originally given by Wainscoat & Cowie 
(1992) in that the peak transmission of our filter is slightly 
redder. Transmission data for these IR filters were mea- 
sured at 77 K. The traces were then multiplied with a 
theoretical model for atmospheric transmission generated 
by the ATRAN program (Lord 1992) as atmospheric ab- 
sorption can be significant at the edges of the bandpasses. 
We also included the (negligible) effect of the quantum 
efficiency (QE) of the IR detectors in these bandpasses. 

The HST bandpasses include the filter transmission pro- 
files, the spectral response of the HST mirror and instru- 
ment optics, and the QE of the detectors. The latter is 
especially important for the NICMOS filters as the QE 
rises sharply from 1-2 i^m (MacKenty et al. 1997) due to 
the low operating temperature of the instrument's HgCdTe 
arrays. 

Also, we compute the ordinary F — 7c, V — K, and J — K 
integrated colors of the SSPs and the strengths of several 
absorption fines (i//3, Mg2, Mgb, H7A, and C4668) de- 
fined on the standard Lick/IDS system (Worthey 1994; 
Worthey et al. 1994; Worthey & Ottaviani 1997). We use 
the analytic fitting functions derived in these studies for 
index strength as a function of stellar temperature, grav- 
ity, and metallicity. We note that a-element enhancement 
in massive elliptical galaxies leads to [Mg/Fe] ratios above 
that of the most metal-rich stars in the solar neighborhood 
(Worthey et al. 1992). Since our models use scaled-solar 
metallicity, we expect our predictions for the Mg2 indices 
to be inaccurate and tabulate them only for completeness. 
Elemental enhancements, however, should have little effect 
on the other spcctrophotometric properties of the model 
stellar populations at fixed total metallicity (A. Bressan 
et al., in preparation). 

3. RESULTS 

3.1. SBF Magnitudes 

We compute M for SSPs as a function of age and metal- 
licity using several combinations of evolutionary tracks 
and SEDs. In rough terms, models with Z > 0.008 and 
ages > 5 Gyr can be thought of as corresponding to the 
classical picture of old, metal-rich stellar populations in 
giant elliptical galaxies, while lower metallicity but com- 



parably old SSPs would correspond to Galactic globular 
clusters. 

Table 2 presents the SBF magnitudes computed using 
the Padova evolutionary tracks and the semi-empirical 
SEDs. Given the wide metallicity range spanned by the 
Padova tracks and the corrections applied by LCB97 to 
make the semi-empirical SEDs agree with observations, 
we consider this combination of tracks and SEDs as our 
"standard models." 

A more enlightening representation of our results is 
given in Figure 3, which shows contour plots of M as 
a function of age and metallicity for a subset of band- 
passes computed with our standard models. For com- 
parison, Figure 4 shows contour plots for a set of inte- 
grated galaxy colors and absorption indices as predicted by 
the same models. Any quantity (SBF magnitude, galaxy 
color, or absorption index) with all vertical contours would 
be completely age-dependent, whereas any quantity with 
all horizontal contours would be completely metallicity- 
dependent. Figure 4 shows that the broad-band colors and 
absorption indices have contours which are nearly straight 
lines. This is the origin of the "3/2 rule" of Worthey 
(1994), which states that changes by a factor of two in 
metallicity roughly mimic changes by a factor of three in 
age for the integrated colors and line strengths of old stellar 
populations. Figure 4 does reveal some differences in the 
slopes of the contours for different integrated colors and 
lino strengths, which reflects the varying age and metallic- 
ity dependences. For example, H/3 is mostly age-sensitive 
while C4668 is mostly metallicity-sensitive. 

Although the contours of the fluctuation magnitudes are 
not as regular as those of the integrated colors and line 
strengths. Figure 3 does reveal regions in age and metal- 
licity where the behavior of M is expected to be relatively 
simple. The optical SBFs of more metal-rich {Z > 0.004) 
populations are largely age-independent, as are near-IR 
(JHK) SBFs for i > 5 Gyr. This is the age/metalficity 
regime expected to be relevant for bright elliptical galax- 
ies. At longer IR wavelengths (L-band and redward), the 
SBFs are mostly agc-clependent. However, we caution that 
these predictions rely more heavily on our prescription for 
stars in the superwind phase on the TP-AGB (which have 
K — L > 1), and the stellar spectral libraries are less ro- 
bust at these wavelengths — our models have not been 
well-tested at wavelengths A > 2.4 /xm due to the lack of 
strong observational constraints. 

3.2. Fractional Contributions of Different Evolutionary 

Phases 

It is instructive to examine the fractional contributions 
to M from different phases of stellar evolution (main se- 
quence, sub-giant branch [SGB], RGB, horizontal branch 
[HB], AGB, bare planetary nebula nucleus, and white 
dwarf). These can be computed from the corresponding 
terms in the numerator of equation (1). Such valuable in- 
formation cannot be currently acquired from observations 
given the challenges of resolving, let alone categorizing, 
individual stars in galaxies. This is especially true for the 
old metal-rich populations of ellipticals because there are 
few nearby examples. Theoretical predictions can indicate 



generically to the SBF absolute magnitude in any filter bandpass. This could also represent the SBF absolute magnitude in the M-band filter 
(Ae// = 4.8 //m), but we avoid using this representation except in Table 2. 
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the sensitivity of M to different stages of stellar evolution 
and, consequently, help to identify which aspects of stellar 
evolution theory can be tested via SBF observations. 

The relative contribution of each evolutionary phase to 
the fluctuation magnitude in a given bandpass is a function 
of both the age and metallicity of the stellar population. 
For f 3 Gyr, the relative contributions of the differ- 
ent phases depend mostly on metallicity, staying roughly 
constant as the population ages. Figure 5 gives a con- 
tour plot of the fractional contribution of AGB stars to 
BVIcJKLasa. function of age and metallicity for our 
standard models. The IR (A > 1 /xm) SBFs arise almost 
exclusively from the RGB and AGB stars, with the AGB 
stars dominating at all mctallicitics for t ^ 3 Gyr. 

Optical SBFs depend more significantly on evolutionary 
phases other than the RGB and AGB. In the Jc-band, the 
HB contribution rises with metallicity, from about 1% at 
Z = 0.0001 (2'0/2OO) to about 20% at Z = 0.05 (2.5 Zq) 
for a > 5 Gyr old population. It is worth noting that in 
our models, AGB stars are important contributors to even 
the /-band SBFs, unlike in the Worthey (1993a) models. 
In the V and i?c-band, the HB contribution is higher still 
at the expense of the RGB and AGB, and in the B-band, 
planetary nebula nuclei become important for Z ^ 0.02 
(Zq). For comparison, a decomposition of the integrated 
light of a solar-metallicity SSP from each evolutionary 
phase is given in Figure 2 of Bruzual & Chariot (1993) 
and in Bruzual (1996). 

3.3. Calibration of SBFs for Distance Determinations 

We now use our models to calibrate M as a distance in- 
dicator via equation (1). Since the ages and metallicities of 
the observed galaxies are not known a priori, the hope is to 
calibrate variations in M arising from variations in these 
parameters by using another observable quantity, such as 
the integrated galaxy color or spectral index. Tonry et al. 
(1997) have derived an accurate calibration for Ic', after 
applying a correction based on a galaxy's V ^ Ic color, 
they find the intrinsic scatter in Ic to be only 0.05 mag. 
Preliminary data for K' reveal no obvious correlation with 
V — Ic color, though the K' sample is still too small for 
a firm conclusion (Jensen et al. 1998). Given the difficul- 
ties involved in tying the SBF observations to the Cepheid 
distance scale (see discussion in Tonry et al. 1999), it is 
important to determine the theoretical zeropoint of the Ic 
versus V — Ic relatiom 

Figure 6 shows Ic J K as a. function of some integrated 
colors and absorption-line indices for models with ages of 
3-17 Gyr and Z = 0.0004, 0.004, 0.008, 0.02, and 0.05. 
Fluctuation magnitudes appear to be the best behaved in 
the Jc-band, i.e., variations in Ic with age and metallic- 
ity are well-correlated with the integrated properties. To 
gauge how different observables are correlated with M in 
different bandpasses, we performed a robust linear fit to 
the fluctuation magnitudes as a function of the various col- 
ors and line indices for the 5-17 Gyr and Z = 0.004 — 0.05 
models. We then calculated the rms scatter of the residu- 
als to these fits. The rms scatter indicates the effectiveness 
of a particular color or index in compensating for stellar 
population variations. The results of these fits are given 
in Table 3. We compare them with observations in § 4.2. 

The scatter in Ic in Figure 6 is only rs 0.1 mag (i.e., 5% 



in distance) as a function oiV — K or J — K color, while it 
is 0.19 mag as a function ofV — Ic- This suggests that the 
accuracy of Ic measurements could be improved by using 
V — K as a^alibrator instead ofV—Ic as is currently done. 
Moreover, Ic as a function of V_— K has a shallower slope 
(1.4, as compared to 4.6 for Ic as a function of V — Ic) 
meaning a greater tolerance of photometric and reddening 
measurement errors. Finally, the model-predicted scatter 
for Ic is also typically smaller than for that for J and K. 

Spectral features can also be used to calibrate SBF 
magnitudes (Thomsen et al. 1997), though uncertainties 
in the effects of a-element enhancement are a concern 
for the models (§ 2.3.2). The metallicity-sensitive in- 
dices Mg2 and C4668 generally have the smallest residuals 
(« 0.1 — 0.2 mag) with 11(3 having the largest residuals. 
This is expected since SBFs of > 5 Gyr populations are 
mostly driven by metallicity. 

Finally, the results of Figure 6 suggest that the SBFs 
measured with a filter around 1 /xm, i.e., between the Ic 
and J bandpasses, should be nearly independent of both 
age and metallicity (leading to horizontal lines in these 
plots). This was also found by Worthey (1993a) who used 
very different models, suggesting the conclusion is robust. 

3.4. k-Corrections 

Optical and near-infrared SBFs can be detected to 

cz ^ 10,000 km s^^ using both groimd-bascd telescopes 
and HST (e.g., Lauer et al. 1998). To derive accurate dis- 
tances to galaxies, and hence an accurate measure of Hq, 
one first needs to apply fc-corrections to the observed SBF 
magnitudes and integrated colors of the galaxies. Such 
corrections can be derived only from models. 

We compute fc-corrections for the integrated galaxy col- 
ors and for M in various bandpasses useful for observ- 
ing distant galaxies. We calculate corrections out to 
cz = 15,000 km s"^ {z = 0.04) for models with Z/Zq = 
{0.2,0.4,1.0,2.5} and ages of 5, 8, 12, and 17 Gyr. The 
corrections for the SBF magnitudes depend mostly on the 
metallicity and very little on the age. They correlate al- 
most linearly with rcdshift over the range considered. Ta- 
ble 4 lists the resulting fc-corrections as a linear function 
of redshift for each metallicity. Wc calculated the slope of 
k{z) for each of the four model ages, and the average and 

rms of these slopes are tabulated. The sign convention 

used for the fc-correction is fc(z) = M{z) — M{z ~ 0), i.e., 
the term is subtracted from observations at z > in order 
to compare to local observations (Humason et al. 1956; 
Rowan- Robinson 1985). 

The amplitudes of the fc-corrections can be significant. 
For instance, at distances of cz « 10,000 km s"'^, the fc- 
correction computed for Ic is ^0.2 mag with little depen- 
dence on the choice of age and metallicity. This corre- 
sponds to a 10% correction in distance (and hence Hq). 
Also, our fc-corrections for Ic and V — Ic agree well with 
those of Tonry et al. (1997). For the IR bandpasses, the 
fc-corrections have a stronger dependence on metallicity, 
but the amplitudes of the corrections are generally smaller 
than in the optical. There are some small irregularities 
seen in the fc-corrections for the 2 /im filters: for K' and 
ii', the fc-corrections can have different signs depending on 
the choice of metallicity, though the amplitudes are small. 
Also, the fc-corrections for the K and K' filters differ sig- 
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nificantly despite their small wavelength separation. This 
is due to the CO bandhead at 2.3 nm, which appears in ab- 
sorption in the spectra of early-type galaxies (Frogel et al. 
1975, 1978). Our /T-band fc-corrections are significantly 
larger than those used by Jensen ct al. (1998), which were 
computed from the Worthey (1994) models, though the ac- 
tual amplitudes of our corrections are still relatively small. 
(Note that the Jensen et al. corrections have the opposite 
sign convention.) 

We recommend using the /c-corrcctions derived from the 
Z = 0.008 and/or Z = 0.02 models, since these models 
roughly span the SBF observational constraints discussed 
below. 

4. COMPARISON WITH OBSERVATIONS AND 
CONSTRAINTS ON STELLAR POPULATIONS 

4.1. SBFs of Globular Clusters 

SBF measurements of globular clusters arc useful tests 
of the models for two main reasons: (1) these systems are 
believed to comprise stellar populations with internally ho- 
mogenous mctallicity and age, and (2) they extend to lower 
metallicities than the available galaxy sample. The only 
observations available are those of Ajhar & Tonry (1994). 
Figure 7 shows their V and Ic data as a function of [Fe/H] 
compared to our standard models.^ We have updated the 
distances used by Ajhar & Tonry (1994) to the Hipparcos 
distance scale from Carretta et al. (2000); the net result 
is to make the SBF magnitudes wO.35 mag brighter than 
reported in the original paper. The agreement between 
models and data is generally good, though the models 
are perhaps too faint by « 0.2 mag for [Fe/H] « -0.7 
{Z = 0.004). At higher metallicities, the models also seem 
to be too red in V-Ic by about the same amount. (See 
also § 4.4.) 

4.2. Ic Measurements of Galaxies 

Since the /-band SBF data comprise by far the largest 
set of measurements to date, they serve as the most_ strin- 
gent test of the models. The calibration of the Ic zero- 
point by Tonry et al. (1997) relied on observations of early- 
type galaxies belonging to seven galaxy groups with HST 
Ccphcid distances, and a much larger sample of galaxies 
was used to determine the slope of the Ic versus V—Ic rela- 
tion. In Figure 8, we show the empirical calibration for Ic 
as a function of V ^ Ic galaxy color derived by Tonry et al. 
(1997). This relation has been re-calibrated by Tonry et al. 
(1999) based on the bulges of six spiral galaxies with both 
SBF and HST Cepheid distances; the numerical results are 
the same as their initial determination. A slightly different 
calibration of the same Ic data was derived by Ferrarese 
et al. (1999), leading to a 0.05 mag brighter zcropoint. 

Also shown in Figure 8 are the results of our standard 
models (Padova tracks with semi-empirical SEDs) and 
those of models with the Pado\^ tracks but with the theo- 
reteical SEDs. The results for Ic from our standard models 

'^The values of [Fe/H] for the data in Figure 7 are those tabulated in Ajhar & Tonry (1994), which are from Zinn (1985). Recently Carretta 
& Gratton (1997, hereinafter CG97) have obtained high-dispersion spectra for a subset of the Zinn (1985) clusters to derive more accurate 
metallicities. They offer a quadratic equation for transforming from the [Fe/H] of Zinn (1985) to their measurements. The effect of changing to 
the CG97 scale is to shghtly increase (<5[Fc/H] » 0.1 - 0.2) the mctallicity of clusters with [Fe/H] < -1.0 and to slightly decrease ((5[Fe/H] « 0.1) 
the metallicity of more metal-rich clusters. The net effect on Figure 7 is small: the [Fe/H] range of the data is slightly compressed, and the 
bulk of the points are shifted by 0.15 — 0.20 dex to higher [Fe/H]. 

*http://199. 120.161. 183/"'worthey/dial/dial_a_model. Note that this version updates the models in Worthey (1994) by including a more 
realistic treatment of the horizontal branch for [Fe/H]< —1.0. 



are in excellent agreement with the empirical calibration. 
The slope fitted to the models (§ 3.3 and Table 3) is 4.56, 
to be compared with the value 4.50 ± 0.25 measured. The 
model-predicted zeropoint at V — Ic = 1.15 is -1.79 mag, 
to be compared with the value of 1.74±0.07 mag de- 
rived empirically by Tonry et al. In contrast, the mod- 
els with theoretical SEDs miss the observational slope by 
« 0.05 — 0.10 mag in V — Ic color. Therefore, it appears 
that the corrections introduced by LCB97 to match the 
observed color-temperature relations of local stars are also 
useful in improving the agreement with Ic observations. 

We note, however, that th e slope predicted from our 
standard models for F814Vt^ as a function oi V — Ic is 
4.42, while the slope measured by Ajhar ct al. (1997) from 
16 galaxies is 6.5 ± 0.7. (See discussion in Ferrarese et al. 
1999 about this point.) A comparison (not plotted) with 
the F814W measurements also shows our standard models 
match better than models using the theoretical SEDs. 

Figure 8 shows thatboth Ic and V — Ic vary with age and 
metallicity, although Ic depends more strongly on metallic- 
ity than age. At the very lowest metallicities (Z < 0.004), 
the model Ic magnitudes flatten out for V — Ic < 1.1. This 
reflects the constant tip magnitude and small color changes 
of the RGB as a function of metallicity at low [Fe/H], 
which is seen in Galactic globular clusters (Da Costa & 
Armandroff 1990). Tonry et al. (1997) observed such a 
deviation in the Ic data for the dwarf elliptical_M 110 
(NGC 205); regions with F - 7^ « 0.8 - 1.0 had Ic fainter 
than expected from extrapolating the linear relation. 

For comparison, we also show in Figure 8 SBF predic- 
tions by Worthey (1994) and Worthey (1999). We used 
Worthey's on-line interpolation engine^ to interpolate his 
results onto our metallicity grid. One striking property 
of these models is the small spread in the model locus, 
implying that variations in age and metallicity are very 
degenerate. The good agreement of the Worthey (1994) 
[Fe/H] = -0.50 to -1-0.50 models with the observations was 
taken by Tonry et al. (1997) as validation of their empir- 
ical calibration. However, it is worth pointing out that 
the Worthey (1994) results for Ic are too faint compared 
to the data for V — Ic ^ 1.1, imhke our models. For 
y — Jc > 1.1, the Worthey (1994) models agree fairly well 
with our models, though there are notable differences in 
the values of Ic and V — Ic from the two set of models for 
any given age and metallicity. Also shown in Figure 8 are 
the Worthey (1999) models, which are based on the same 
Padova tracks as our standard models for stars from the 
main-sequence through the early AGB. The results from 
his updated models differ significantly from our own and 
those of Worthey (1994). Moreover, they do not match 
the data very well. 

4.3. K' Measurements of Galaxies 

In Figure 9, we compare our standard models to K' data 
from the literature. The ranges of ages and metallicities in 
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the models are the same as in Figure 8. The observations 
comprise the good S/N measurements of Fornax, Virgo, 
and Eridanus galaxies from Jensen et al. (1998), using up- 
dated 7-band SBF distances to the galaxies from Tonry & 
Blakeslcc (2000); this sample is composed almost entirely 
of luminous galaxies {Mb < —20). We do not include the 
lower S/N data of Jensen et al. given the concerns about 
systematic biases in SBF measurements at low S/N (see 
Jensen et al. 1996). For the same reason, we do not include 
data from Pahre & Mould (1994), and also because these 
authors did not correct their SBF measurements for glob- 
ular cluster contamination. Finally, we include K' data 
for M 31 and M 32 from Luppino & Tonry (1993). 

Unlike the Ic data in Figure 8, Figure 9 shows that the 
K' data do not strongly favor our standard models over 
the models with the theoretical SEDs. Also in contrast to 
Ic, the effects of age and metallicity in the {V — Ic, K'} 
plane are more distinct, as can be seen by the roughly or- 
thogonal orientations between models with common ages 
and those with common metallicities. (This effect was also 
noted by Worthey 1993a). The K' models become slightly 
fainter with increasing age. This is expected from the core 
mass-luminosity relation for AGB stars. 

The observations in Figure 9 span a range of model ages 
of about a factor of three, and the galaxies with the bright- 
est if'-band SBFs have SSP-equivalent ages near 3 Gyr. 
The inferred metallicities arc roughly solar, which is some- 
what higher than the value of w 0.5Z© inferred by Jensen 
et al. (1998) using the Worthey (1994) models.^ At these 
ages and metallicities, the brightest AGB stars are M- 
type stars. The data in Figure 9 also show that the bluer 
( V — Jc < 1.2) galaxies in the sample tend to have brighter 
K' . The implication from the models is that the stellar 
populations dominating the K'-hand SBFs in the bluer 
galaxies arc younger than those in the redder galaxies; in 
addition, the bright K' of the bluer galaxies suggests they 
also have roughly solar metallicities. 

4.4. Optical/IR Fluctuation Colors of Galaxies and 

Globular Clusters 

SBF color measurements provide another interesting 
test of the models. Since the fluctuations are almost en- 
tirely dominated by RGB and AGB stars, SBF colors are 
a better tracer of the giant branch colors than the inte- 
grated light. In particular, imlike the integrated colors 
used above, SBF colors are independent of model uncer- 
tainties in the main sequence, subgiant branch, and hor- 
izontal branch stars. Also from a practical standpoint, 
SBF colors are distance-independent (neglecting the k- 
corrcctions). Thus they avoid any systematic errors from 
uncertainties in the galaxy distances, which are inherent 
in using absolute fluctuation magnitudes. However, the 
number of galaxies to date with multi-color SBF data is 
relatively small. 

We compare our standard models with galaxies with 
good optical/IR SBF color data in Figure 10. The data are 
for M 31, M 32, and Virgo cluster galaxies. The IR data 
are the same as in Figure 9, and the optical data arc from 
Tonry et al. (1990), with some minor revisions from Tonry 

^Kuntschner (2000) used single-burst population synthesis models 
galaxies in Figure 9 (NGC 1379, 1399, and 1404). The ages and n 
comparable to his findings. 



& Blakeslee (2000). For the V Rc Ic colors, most of the 
models overlap each other (Figure 10a). This happens be- 
cause of the turnover of the tip of the giant branch and the 
increasing contribution from core-He burning stars with in- 
creasing metallicity (Worthey 1993a; Ajhar & Tonry 1994). 
Therefore, these colors are largely _degener ate to age and 
metallicity. The optical/IR {V Ic K') model colors (Fig- 
ure 10b) show the same reversal as described above in 
the V-Ic color, but the Ic-K' color increases monotoni- 
cally with metallicity. This plot clearly demonstrates that 
the Ic^K' color is predicted to be very metallicity-sensitive 
and highly age-insensitive. We return to this point in § 5.3. 

There is generally good agreement between our standard 
models and the data. The comparison with observed opti- 
cal SBF colors in Figure 10a indicates that our predictions 
for V-Ic are consistent with the data, while the predicted 
Rc-Ic color may be «0.1 mag too blue. However, Fig- 
ure 10b would instead suggest that the model V-Ic colors 
are too red by «0.15 mag, an offset which is also sug- 
gested by comparisons with globular cluster SBF colors 
(§ 4.1). Thus, the models have some difficulty in match- 
ing all the observed SBF colors simultaneously to better 
than 0.1 0.15 mag. Blakeslee et al. (1999) did a similar 
comparison using the Worthey (1994) models and found 
a disagreement of about 0.3 mag between the V-Ic model 
colors and the data, which they suggested was evidence 
for multi-metallicity stellar populations in ellipticals. Our 
models do not show this large a discrepancy. Furthermore, 
we point oiit that the relatively small wavelength lever- 
age of the V-Ic color is not a very strong discriminant 
— the B-Ic color should more sensitive to populations 
of composite metallicity in ellipticals. Past constraints on 
such populations from V Rc Ic colors alone (Buzzoni 1993; 
Worthey 1993a) could likely be improved by expanding the 
wavelength range considered. 

The agreement between observations and models is 
worsened if we use theoretical SEDs in our models instead 
of semi-empirical ones. We find that for Rc-Ic and V- 
Ic, the theoretical SEDs produce model SBF colors «0.2 - 
0.3 mag bluer relative to the semi-empirical SEDs. In ad- 
dition, the Ic-K' colors are wO.2 mag too blue compared 
to the data. 

Even though the number of available measurements is 
small, the SBF colors have relatively little scatter from 
galaxy to galaxy. The small observed ranges of Ic-K' 
and V-Ic colors (Figure 10b) suggests that the aperture- 
averaged metallicities of the brightest giant stars are 
roughly consistent from galaxy to galaxy. Also, the Local 
Group galaxies (M 31 and M 32) have SBF colors compara- 
ble to those of the Virgo galaxies, despite the very different 
physical apertures used in the measurements, since Virgo 
is about 20 times farther away than M 31 and M 32. This 
implies that the effects of any SBF color gradients in these 
galaxies are small. 

In Figure 11, we compare the V-Ic and V — Ic colors 
from our standard models with observations of both galax- 
ies and globular clusters. The globidar cluster data are the 
same as in § 4.1 and Figure 7. The galaxy data are from 
Ajhar & Tonry (1994) with some revisions from Tonry & 

interpret the optical absorption-line strengths of three of the Fornax 
allicities of these galaxies inferred in Figure 9 from our models are 
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Blakcsloc (2000); the galaxies have Mb « -18.5 to -21.5 
and belong mostly to the Virgo cluster. The reddening 
corrections adopted for the galaxies have been changed to 
those of Schlegel et al. (1998). Overall, Figure 11 shows 
good agreement between our standard models and the data 
and illustrates the expected metallicity dichotomy between 
Galactic globular clusters and elliptical galaxies. Interest- 
ingly, the bulk of the galaxy data suggests a large range in 
age with a small range in metallicity (mostly sub-solar), 
which is reminiscent of the trends found in the smaller K' 
dataset. For the globular clusters, the V — Ic measurement 
errors are much larger than for the galaxies, and the ob- 
served color spread is consistent with little age dispersion. 

5. DISCUSSION 
5.1. Systematic Effects of Model Uncertainties 

In trying to extract information about the stellar pop- 
ulations of galaxies using population synthesis models, it 
is important to explore the cfFccts of changing the model 
inputs. We have shown that our standard models can 
reproduce most of the current observations, but some 
slight offsets between the models and data (e.g., Figs. 7 
and 10) could be evidence either for systematic errors in 
the model inputs or for real physical differences between 
actual galaxy populations and SSPs. 

The freedom of choice in the evolutionary tracks 
(Padova or Geneva) and spectral libraries (empirical, semi- 
empirical and theoretical) of our models allows us to ex- 
plore the range of uncertainties in our SBF predictions. 
There are appreciable differences in the predicted fluctu- 
ation magnitudes depending on the choice of model in- 
puts, though the differences do not follow obvious trends 
in metallicity, age or bandpass. 

5.1.1. Multi- Metallicity Models: Choice of SEDs 

In comparing the results of multi-metallicity models 
from the Padova tracks combined with either the semi- 
empirical or theoretical SEDs, the most noticeable differ- 
ences arise in the /-band region [Ic and FSIAW). For 
all metallicities Z > 0.004, Ic obtained from the semi- 
empirical SEDs are about 0.1-0.2 mag brighter than those 
obtained from the theoretical SEDs. Similar differences 
occur over smaller metallicity ranges for J {Z > 0.02) and 
H {Z = 0.004 — 0.008), but with the semi-empirical predic- 
tions now being fainter than the ones from the theoretical 
SEDs. Since the predicted SBF colors therefore differ by 
0.2-0.4 mag in this metallicity range, J and Ic-H data 
could provide an independent test of the LCB97 SED cor- 
rections (§ 2.3). The SBF magnitudes from the two sets 
of SEDs otherwise agree at the <0.1 mag level. Similarly, 
the integrated colors and line indices obtained with both 
sets are comparable, with the semi-empirical SEDs giving 
slightly redder V — Ic colors than the theoretical ones. 

5.1.2. Solar-Metallicity Models: Choice of Evolutionary 
Tracks and SEDs 

Table 5 provides the SBF magnitudes, integrated col- 
ors, and line indices from the six possible types of solar- 
metallicity models we can investigate from our two sets of 
evolutionary tracks (Padova and Geneva) and three sets of 
SEDs (theoretical, semi-empirical, and empirical). Note 



that, by construction, the theoretical SEDs should pro- 
duce less accurate results than the semi- empirical SEDs, 
which have been adjusted to match observations of solar- 
neighborhood stars. There are considerable (> 0.3 mag) 
variations in the fluctuation magnitudes from different 
choices of tracks and spectral libraries. For a given set of 
evolutionary tracks (Padova or Geneva), the predictions of 
the three different sets of SEDs are relatively consistent, 
with differences at the <0.15 mag level in all bandpasses 
(except for the Af-band SBFs from the empirical SEDs). 

However, for a fixed set of SEDs, larger differences arise 
when changing the stellar evolution prescription. The 
fluctuation magnitudes from the Geneva models are up 
to 0.3 mag brighter in the optical than those from the 
Padova models and w 0.2 — 0.5 mag fainter in the in- 
frared. This is caused by the significant difference in 
the relative numbers of RGB and AGB stars predicted 
by the two tracks. (See also discussions in Bruzual 1996 
and Chariot et al. 1996.) The Geneva models have a 
much larger AGB contribution to the SBF signal at solar- 
metallicity, especially from the early AGB phases. The 
result is that the Geneva tracks produce optical SBFs 
somewhat brighter and near-IR SBFs somewhat fainter 
than the Padova tracks. In contrast to these large differ- 
ences in SBF predictions, the integrated properties from 
the Padova and Geneva models at solar metallicity are 
very comparable. Bruzual et al. (1997) examined the in- 
tegrated spectra of near solar-metallicity Galactic Bulge 
globular clusters and found good agreement using semi- 
empirical SEDs with either the Padova or Geneva tracks. 

Figure 12 compares all six solar-metallicity models 
against the data. In general, for a fixed set of stellar 
evolutionary tracks, the theoretical SEDs produce bluer 
integrated and SBF colors as well as a slightly fainter Ic; 
the semi-empirical and empirical SEDs give very compa- 
rable results, which is reassuring. For a fixed choice of 
input SEDs, the Geneva tracks produce slightly redder in- 
tegrated and optical SBF colors than_the Padova tracks. 
The Geneva tracks also give brighter Ic and fainter K and 
hence bluer Ic-K. 

Overall, the SBF magnitudes (Figure 12a and b) fa- 
vor the Padova models, especially the Ic data. The con- 
straints from SBF colors are less conclusive. For the op- 
tical colors {V Rc Ic), the choice of evolutionary tracks 
and SEDs are largely degenerate (Figure 12c), with the 
solar-metallicity models occupying a locus similar to that 
of the multi-metallicity models in Figure 10. For the 
optical/near-IR colors {V Ic K), the slight discrepancy 
noted in § 4.4 between our standard models and the data 
is not found when adopting the Geneva tracks with the 
theoretical SEDs (Figure 12d), though this combination is 
largely disfavored by the other comparisons in Figure 12. 

5.2. Constraints on the IMF and Stellar Evolutionary 

Lifetimes 

As noted in_§ 4.2, the tightness of the observed correla- 
tion between Ic and V — Ic for elliptical galaxies reflects 

a degeneracy between the effects of age and metallicity on 
these observables. This degeneracy may be useful to con- 
strain variations in other parameters of the models, such 
as the IMF or the uncertainties in the evolution of the cool 
stars. 
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The tightness of the observed relation between Ic and 
V — Ic sets a limit on how different the IMF of elliptieal 
galaxies can be from that in the solar neighborhood. Fig- 
ure 13 shows the effect of changing the IMF slope on the 
correlation between /<, and V — Ic- Lowering the IMF slope 
relative to Salpeter has only a small effect; therefore the 
data are not very sensitive to IMFs flatter than Salpeter. 
A similar conclusion has been reached when purely us- 
ing integrated colors to constrain the IMF (e.g. Tinsley 
1978; Frogel 1988). However, Figure 13 does suggests that 
the IMF in observed elliptical galaxies cannot be much 
steeper (i.e., more dwarf-dominated) than the IMF in the 
solar neighborhood. This interesting constraint from SBF 
data complements those on the higher-mass end of the 
IMF from chemical abundance and integrated light stud- 
ies (e.g., Worthey et al. 1992; Kennicutt 1998). 

Since the evolutionary lifetimes of low- and 
intermediate-mass stars are still a major uncertainty in 
population synthesis models (Chariot et al. 1996), it is 
also revealing to investigate how these are constrained by 
the tight correlation between Ic and V—Ic color. Figure 14 
illustrates the effect of changing the lifetimes of post-main 
sequence evolutionary phases (RGB, HB, and AGB) by a 
global factor in our standard models. The principal effect 
is to change Ic, with small changes in the V — Ic color. 
This comparison indicates that model uncertainties in the 
lifetimes of the post-MS evolutionary phases are probably 
less than ±50%. (Conversely, any large systematic change 
to the Cepheid distances which set the zeropoint for the 
/-band SBF calibration would suggest that the model life- 
times need to be adjusted.) This constraint is interesting 
because it applies at the relatively high metallicities of 
elliptical galaxies. For low-metallicity stellar populations, 
useful constraints on the accuracy of the post-MS life- 
times can be derived from star counts in Galactic globular 
clusters (e.g., Renzini & Fusi Pecci 1988). 

5.3. New Tools for Breaking the Age/Metallicity 
Degeneracy 

Broad-band colors are largely degenerate to changes in 
age and metallicity in old populations, with changes of 
S = d{\og age)/d{\og Z) « 3/2 approximately preserving 
the colors (Worthey 1994). Stellar absorption-line indices 
can be more sensitive to either age or metallicity, e.g., H(i 
and are age-sensitive {S < 1.0) and Mg2 and C4668 
are metallieity-sensitive (S* w 2 — 5). SBFs are expected 
to depend mostly on metallicity in old populations be- 
cause they closely track the temperature of the RGB and 
AGB, whose colors are governed by metallicity (§ 2.2). 
SBF magnitudes and colors of old populations at near- 
IR wavelengths are predicted by our models and those of 
Worthey (1994) to have d(log a5e)/rf(log Z) > 5 - 6. This 
strong sensitivity to metalhcity suggests that SBF data, 
when combined with age-sensitive observables, could ef- 
fectively disentangle the effects of age and metallicity in 
interpreting unresolved stellar populations. 

A full investigation of using SBFs for stellar population 
studies is beyond the scope of this paper; in addition, the 
existing datasets for such analyses are limited. Instead we 
highlight in Figure 15 two possible methods to break the 
age-metallicity degeneracy, both of which rely on the Ic-K 
color as a metallicity indicator. Since Balmer absorption 



lines are standard age indicators, the use of II/3 in com- 
bination with Ic-K should be effective in distinguishing 
age from metallicity (Figure 15a). Similar results are ob- 
tained when using Il7^ as an age indicator. Our models 
also predict that L (3.5 /im) and M-band (4.8 /im) SBFs 
are very sensitive to age (Figure 15b), though this predic- 
tion should be treated with caution (see § 3.1). Given the 
very high thermal background of the Earth's atmosphere, 
L-band SBFs are measureable from the ground only for 
the very closest galaxies su(;h as M 31, while observing 
M-band SBFs is presently impossible. The future Space 
Infrared Telescope Facility (SIRTF) will enable SBF mea- 
surements in both of these bands. 

SBF colors such as Ic^K could offer some advantages 
over metal-absorption lines as metallicity indicators for 
elliptical galaxies. The interpretation of metal lines is 
generally hampered by some possible complications: (1) 
the effect of selective a-element enhancement in elliptical 
galaxies which directly affects indices such as Mg2, and 
(2) the limited range in stellar temperatures and gravities 
over which the analytic fitting functions have been derived 
to parameterize line strengths (e.g., Worthey et al. 1994). 
Also, absorption lines are typically measured only in the 
central regions of galaxies, while SBF measurements sam- 
ple a much larger area. Finally, as seen in Figure 15, be- 
cause of their weighting to the most luminous cool stars, 
SBFs offer a much greater dynamic range than integrated 
spectral properties; changes in metallicity should be more 
clearly detectable in SBF measurements. 

6. CONCLUSIONS 

We have presented theoretical predictions for SBFs of 
single-burst stellar populations (SSPs) spanning a wide 
range of ages (from 1 to 17 Gyr) and metallicities (from 
1/200 to 2.5 times solar). Our calculations are based 
on the population synthesis models of Bruzual & Char- 
lot (2000), in which the stellar evolution prescription and 
spectral libraries are improved over the models used in 
previous SBF studies. In particular, our models have 
been optimized during the course of this work by refin- 
ing the prescription for the latest phases of stellar evo- 
lution, which are important contributors to the optical 
and infrared SBF signal. Our standard predictions are 
based on multi-metallicity evolutionary tracks from the 
Padova school and semi-empirical stellar spectra designed 
to match the observed color-temperature relations of solar- 
neighborhood stars at solar metallicity (LCB97). 

Using our models, we generate several basic predic- 
tions as a function of age and metallicity. 

1. We compute SBF magnitudes and integrated colors 
for a large set of ground-based and space-based 
(HST) optical and infrared bandpasses. These are 
supplemented with the strengths of several optical 
absorption- line indices on the Lick/IDS system. 

2. We provide results for solar- metallicity models 

using several combinations of stellar evolutionary 
tracks and spectral libraries. These can be used 
to assess the systematic effects of different model 
inputs on the results. 

3. We predict the fractional contribution of different 
stellar evolutionary phases to the SBFs. Since 



12 



Liu, Chariot & Graham 



this information cannot be easily derived from 
observations, the models provide insight into which 
phases are important contributors to the SBF 
signal for a given bandpass. 

Our model results directly benefit SBF distance deter- 
minations, specifically: 

4. We use the models to determine purely theoretical 

calibrations for SBFs in many bandpasscs. These 
are independent of any systematic errors in 
Cepheid distances or reddening corrections, which 
affect only empirical calibrations. 

5. Wc tabulate fc-corrections out to 

cz < 15,000 km s^^ (z = 0.04), which are 
required for accurate determinations of Hq. Wc 
find that the A;-corrections are roughly linear in 
this redshift range. Metallicity has a stronger 
effect on fc-corrcctions in the near-infrared than 
in optical, but the amplitudes of the corrections 
are also generally smaller in the near-infrared. We 
conclude that systematic errors from uncertainties 
in the fc-corrections are not important sources of 
error for Hq determinations. 

6. We suggest that the scatter in /-band SBF 
distances can be further reduced by using the 
integrated V — K galaxy color instead oiV — Ic to 
correct for stellar population variations between 
galaxies. The reason for this improvement is that 
the V — K color is more sensitive to metallicity, 
which also drives the Ic signal. 

7. Our models predict that the fluctuation magni- 
tudes should be independent of population age and 

metallicity around 1 /im. A similar conclusion was 
reached by Worthey (1993a) using very different 
population synthesis models — this suggests that 
the prediction is robust. Therefore, observations 
taken with a .Z^-band filter from a large (8-10 m) 
ground-based telescope or with the _F110VF filter 
in HST NICMOS should allow SBF distance 
measurements which are more robust against 
galaxian population variations. 

Wc have compared onv model results with nearly all the 
SBF measurements available to date. Since SBFs are es- 
pecially sensitive to the cool, luminous stars on the upper 
RGB and AGB, they provide important tests for popu- 
lation synthesis models. The existing dataset comprises 
Galactic globular clusters, M 31, M 32, and early-type 
galaxies in nearby clusters. The /-band dataset is by far 
the most extensive; there are some /C-band SBF magni- 
tudes and optical/IR SBF colors, but more measurements 
are needed for further testing. 

We find generally good agreement between mod- 
els and data and also suggest some new tests for the 
models. Specifically: 

8. Our models reproduce V and Ic observations of 
Galactic globular clusters. This test is comple- 
mentary to those based on galaxy data, since the 
globular clusters have much lower metallicities. 



Models with [Fe/H] « -0.7 might be ^0.2 mag too 
red in V-Ic, although more data are needed to 
verify this. 

9. Our standard models provide the best agreement 
to date with the tight empirical calibration of Ic 
over the entire observed range oiV — Ic galaxy 
color. The zeropoint and slope of the calibration 
predicted by our models agree remarkably well 
with those derived from the data. Moreover, the 
models indicate a saturation of /c for V — Ic < 1.0, 
which is also seen in the observations. The reason 
for this flattening is most likely the constancy of 
the /-band tip of the RGB for metal-poor stellar 
populations. The small scatter in the empirical 
calibration as a functicm of V — Ic galaxy color is 
also reproduced by the models; this arises because 
of the partial age/metallicity degeneracy in the 
{Ic, V — Ic} parameter space. This degeneracy is 
a boon for distance measurements. 

10. Our standard models also agree with K' observa- 
tions, although this is based_on a much smaller 
sample of galaxies. In the {/T', V — Ic} parameter 
space, changes in age and metallicity are roughly 
orthogonal. K' brightens for populations of higher 
metallicities and younger ages, as expected from 
observations of the RGB and AGB of star clusters 
in the Galaxy and the Magellanic Clouds. 

11. The optical/IR fluctuation colors predicted by 
our models agree with the observations, although 
some discrepancies exist at the w 0.1 — 0.2 mag 
level. An advantage of testing the models against 
measurements of SBF colors is that the data are 
immune to errors in the galaxy distances. 

12. The semi-empirical SEDs of LCB97 provide 
better agreement with SBF observations at 
all metallicities than their theoretical SEDs. 
Observations of Ic-J and Ic-H colors would help 
to verify this result. For solar metallicity, the 
results obtained from the empirical spectral library 
of Pickles (1998) agree closely with those from the 
semi-empirical library of LCB97. 

13. For solar metallicity, the Padova evolutionary 
tracks seem to provide better agreement with 
SBF observations than the Geneva tracks. The 
integrated spectral properties from the two sets of 
tracks are very comparable at solar metallicities. 
However, the SBF data are a sensitive test 

for deciding between the Geneva and Padova 
tracks, since the differences are larger in the 
SBF predictions than in those for the integrated 
spectra. 

14. From the tightness of the empirical Ic calibration, 
we conclude that the lifetimes of post-main 
sequence phases (RGB, core-He burning, and 
AGB) in the evolutionary tracks are probably 
accurate to within better than ±50%. 
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By comparing our single-burst models with the available 
dataset, mostly composed of luminous galaxies in nearby 
clusters, our preliminary findings on the stellar popula- 
tions dominating the SBFs are: 

15. The metallicities inferred from SBF magnitudes 
and SBF colors show little spread. The metallicities 

favored by the optical/SBF colors arc slightly 
sub-solar, while those favored by the K' data are 
around solar. 

16. SBF color measurements show no obvious 
differences for galaxies obscirvcd with different 
linear aperture sizes, though the available dataset 
is small. The implication is that SBF distance 

measurements should be relatively insensitive to 
systematic errors due to aperture effects. 

17. The _ages_ inferred from comparisons of both K' 
and V-Ic with the V — Ic integrated color span a 
range of about a factor of three, with the youngest 
ones near 3 Gyr. Note that estimates based on 
combinations of SBF colors and integrated colors 
are independent of the galaxy distances. 

18. For old populations, the tightness of the empirical 
/-band SBF calibration also indicates that the 
IMF in elliptical galaxies cannot be significantly 
steeper than that in the solar neighborhood. 

Finally, we suggest that SBF measurements can offer use- 
ful new tools for stellar population studies: 



The single-burst models we have investigated can ac- 
count for the full observed ranges of SBF magnitudes, SBF 
colors, and integrated colors for bright elliptical galaxies 
in nearby clusters. It is important to realize that, al- 
though the SBF observations can be most simply repro- 
duced by models with around solar metallicity and a sig- 
nificant spread in age, a more refined analysis is required 
to interpret these measurements in terms of the star for- 
mation history of elliptical galaxies. In particular, there 
arc multiple lines of evidence that both cluster and field 
elliptical galaxies have experienced more than one episode 
of star formation (e.g., Schweizer & Seitzer 1992; Barger 
ct al. 1996; Poggianti et al. 1999). To constrain the ages 
and metallicities of different stellar generations in elliptical 
galaxies, we then require a combination of various age and 
metallicity indicators. Unfortunately, the published SBF 
and absorption-line studies contain few galaxies in com- 
mon. In a future paper (M. Liu et al., in preparation), we 
exploit a more extensive set of new SBF measurements to 
investigate the stellar content of elliptical galaxies. 
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19. In old populations, the SBF magnitudes and colors 
are predicted to be very sensitive to metallicity, 
especially at near-IR {JHK) wavelengths. 

This may offer a potent means of breaking the 
age/metallicity degeneracy inherent in studies 
based on integrated spectral properties. 

20. We find that the Ic^K SBF color is very sensitive to 
metallicity because of the decreasing temperature 
of the giant branch with increasing metallicity. 
Thus, Ic-K might be used in combination with 
age-sensitive observables such as Balmer absorption 
lines to constrain the ages and metallicities of 
elliptical galaxies. SBF colors may also present 
advantages over metal absorption lines such as Mg2 
and C4668, which are affected by uncertainties in 
the patterns of a-element enhancement in elliptical 
galaxies. 

21. Our models suggest that the L-band and M-band 
SBFs are very sensitive to age, although our 
l)redictions are not optimized in this wavelength 
range. This potentially interesting result should 
be further investigated using more appropriate 
models. 

22. Observations of B-I,, with Ic^K may be usehil to 
identify stellar populations of different metallicities 
in elliptical galaxies. 
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Fig. 1. — Cumulative distribution as a function of evolutionary phase for the integrated light and SBFs of a 12 Gyr old 
solar- metallicity single-burst stellar population from our standard models (Padova tracks with semi-empirical SEDs). 

Evolutionary phases range from the zcro-agc main sequence on the left to the end of the AGB on the right using an 
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Fig. 2. — Theoretical isochrones from the Padova models illustrating the effect of age and metallicity variations for old 
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Fig. 3. — Contours of SBF magnitudes for the BVIcJKL bandpasses in the parameter space of age and metalhcity from 
our standard models (Padova tracks with semi-empirical SEDs). Horizontal contours would indicate that a quantity 
depends only on metallicity, and vertical contours that it depends only on age. 
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Fig. 4. — Contours of predicted galaxy colors and absorption line indices in the parameter space of age and metallicity 
for the same models as in Figure 3. Horizontal contours would indicate that a quantity depends only on metallicity, and 
vertical contours that it depends only on age. 
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Fig. 5. — Fractional contribution of all AGB stars (including carbon stars) to the SBF magnitudes in the BVI^JKL 
bandpasses as a function of age and metallicity for our standard models (Padova tracks with semi-empirical SEDs). For 
the IR (A > 1 \xm) bandpasses, the remaining contribution is almost exclusively from RGB stars. 
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Fig. 6. — Results for UK SBF magnitudes as a function of integrated light properties from our standard models. Ag 
of 3, 5, 8, 12, 17 Gyr are plotted for a range of metallicities, as indicated. Models of a given metallicity have the san 
symbol, with increasing symbol size representing increasing age. 
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Fig. 8. — Comparison of different models with the empirical calibration of Ic as a function oi V — Ic color (heavy line; 
Tonry et al. 1997, 1999). The hatched region represents the la spread in the empirical calibration. Models of a given 
metallicity have the same symbol, with increasing symbol size representing increasing age. Lines connect models with 
the same age. Left: Our model predictions using Padova tracks with semi-empirical (top) and theoretical (bottom) 
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Fig. 9. — Comparison of different models with the K' observations. The models are computed using the Padova 
evolutionary tracks with the semi-empirical SEDs (our standard models; top) and the theoretical SEDs {bottom). The 
high-S/N observations for ellipticals in Virgo, Eridanus, and Fornax (Jensen et al. 1998) are plotted as circles (•), and 
those for M 32 and the bulge of M 31 (Luppino & Tonry 1993) are plotted as stars (*). The colored model symbols are the 
same as in Figure 8. Models with the same metallicity have the same symbol, with increasing symbol size representing 
increasing ages. Lines connect models of the same age at 3, 5, 8, 12, 17 Gyr, with the oldest models having the reddest 
V — Ic colors. 
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Fig. 11. — Comparison of om' standard models with observations of V-Ic fluctuation color as a function of integrated 
V~Ic color. Milky Way globular clusters are shown as squares (□) and nearby early-type galaxies as circles (•). Different 
regions of NGC 205 are shown as stars (★) . For clarity, errors for the globular cluster V ~ Ic colors are shown for only 
one cluster, chosen to have the median error of the sample. Most of the galaxy sample is from the Virgo cluster. 
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Fig. 12. — Comparison of our solar-metallicity models using different combinations of stellar evolutionary tracks (Padova 
versus Geneva) and spectral libraries (semi-empirical [SE], empirical [E], and theoretical [T]) with the observations from 
Figures 8, 9, and 10. All models are shown at ages of 3, 5, 8, 12, 17 Gyr, with increasing symbol size representing 
increasing age. As expected, for a fixed set of evolutionary tracks, the semi-empirical SEDs give results similar to those 
from the empirical SEDs and provide better agreement with the data than the theoretical SEDs (see § 2.3). 
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Fig. 13. — Effect of changing the slope of the IMF in our standard models. The IMF is parameterized as dN — 
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spread in the empirical calibration. The ages and metallicities of the models are the same as in Figure 8. The /-band 
SBF data disfavor an IMF much steeper than Salpeter. 
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Fig. 14. — Effect of changing the lifetimes of the post-main sequence evolutionary phases in our standard models. The 
heavy black line and hatched region represent the la spread in the empirical calibration. The ages and metallicities 
of the models are the same as in Figure 8. The /-band SBF data indicate that the post-main sequence evolutionary 
lifetimes in the models cannot be in error by more than ±50% over wide ranges in age and metallicity. 
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Fig. 15. — Two examples of using SBF measurements to disentangle age and metallicitiy effects in stellar populations 
studies. Top: The age-sensitive H/? absorption line versus the metallicity-sensitive SBF Ic-K color. Bottom: The 
age-sensitive i-band (3.5 /im) SBF magnitude versus the SBF Ic^K color. 
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Table 1 
Model Options 



Inputs 


Options 


Metallicity Range*^ 


Reference 


Evolutionary tracks 


Geneva 






Schaller et al. (1992)'' 




Padova 


{1/200 


- 2.5} X Zq 


Bressan et al. (1993)^ 


Spectral libraries 


theoretical 


{1/200 


- 2.5} X Zq 


Lejeune et al. (1997, 1998) 




semi-empirical 


{1/200 


- 2.5} X Zq 


Lejeune et al. (1997, 1998) 




empirical 




Zq 


Pickles (1998) 



Note. — Our standard models use the Padova tracks with the semi-empirical SEDs. 
^Solar metallicity is Zq = 0.02 in these models. 

•^More complete lists of references are available in Chariot et al. (1996), Bruzual et al. (1997), and Bruzual & Chariot 
(2000). 
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-5, 


.46 


-5, 


,45 


-5, 


.58 


-5, 


,80 


-5, 


,98 


-6, 


.26 


1 


.33 


3 


.34 


0, 


,93 


1, 


,37 


0, 


,29 


4, 


,24 


-7.44 


6.14 



0.0500 


1 


2.26 


1.00 


-0.32 


-2.10 


-1.98 


-4.67 


-4.43 


-5.17 


-6.33 


-6.55 


-7.14 


-7, 


19 


-7.21 


-7. 


,29 


-7.68 


-7.84 


-8.11 





,96 


2.91 


0.98 


3.29 


0.18 


2.53 


1.88 


5.68 


0.0500 


2 


2.86 


1.32 


0.23 


-1.28 


-1.21 


-4.25 


-3.97 


-4.77 


-5.82 


-6.07 


-6.63 


-6, 


,72 


-6.74 


-6 


,91 


-7.30 


-7.53 


-7.86 


1 


,14 


3.22 


0.98 


2.45 


0.23 


3.34 


-3.54 


7.01 


0.0500 


3 


3.02 


1.48 


0.46 


-0.94 


-0.88 


-3.94 


-3.66 


-4.43 


-5.44 


-5.68 


-6.22 


-6, 


,31 


-6.31 


-6, 


.49 


-6.87 


-7.09 


-7.42 


1 


.23 


3.35 


0.98 


2.06 


0.27 


3.83 


-5.41 


7.85 


0.0500 


5 


2.78 


1.67 


0.61 


-0.83 


-0.77 


-3.78 


-3.51 


-4.30 


-5.29 


-5.54 


-6.05 


-6, 


,15 


-6.16 


-6 


.39 


-6.56 


-6.89 


-7.21 


1 


.31 


3.51 


1.00 


1.73 


0.30 


4.35 


-6.91 


8.58 


0.0500 


8 


2.61 


1.77 


0.60 


-0.84 


-0.77 


-3.68 


-3.41 


-4.21 


-5.19 


-5.45 


-5.96 


-6, 


,08 


-6.08 


-6 


.34 


-6.40 


-6.80 


-7.12 


1 


.35 


3.57 


1.00 


1.55 


0.32 


4.66 


-7.64 


8.97 


0.0500 


12 


2.39 


1.75 


0.66 


-0.73 


-0.67 


-3.51 


-3.25 


-4.03 


-5.01 


-5.26 


-5.76 


-5, 


,88 


-5.88 


-6 


.12 


-6.18 


-6.56 


-6.89 


1 


.41 


3.67 


1.01 


1.34 


0.36 


5.02 


-8.49 


9.72 


0.0500 


17 


2.35 


1.82 


0.79 


-0.54 


-0.48 


-3.31 


-3.04 


-3.82 


-4.80 


-5.06 


-5.55 


-5, 


,66 


-5.66 


-5, 


.90 


-5.95 


-6.32 


-6.66 


1 


.44 


3.71 


1.00 


1.19 


0.38 


5.28 


-9.00 


10.33 



*Our standard models use the Padova evolutionary tracks, the semi-empirical SEDs, and a Salpeter IMF. 



Table 3 

Coefficients of Robust Linear Fits to SBF Magnitudes 









Ic 






F81AW 






J 






FWOW 






K 




CO 


Cl 


(T 


CO 


Cl 


a 


CO 


Cl 




CO 


Cl 




CO 


Cl 


a 


as 


a function of V — Ic 


-7.03 


4.56 


0.19 


-6.77 


4.42 


0.18 


-1.95 


-1.45 


0.22 


-2.96 


-1.44 


0.23 


-2.54 


-2.34 


0.31 


as 


a function of V — K 


-5.59 


1.36 


0.10 


-5.38 


1.32 


0.10 


-2.14 


-0.52 


0.17 


-3.12 


-0.53 


0.19 


-2.90 


-0.83 


0.24 


as 


a function of J — K 


-6.50 


5.80 


0.12 


-6.26 


5.62 


0.11 


-1.63 


-2.41 


0.15 


-2.58 


-2.46 


0.16 


-2.12 


-3.78 


0.20 


as 


a function of H/3 


1.05 


-1.41 


0.33 


1.07 


-1.37 


0.32 


-4.29 


0.31 


0.26 


-5.26 


0.30 


0.27 


-6.34 


0.53 


0.38 


as 


a function of Mg2 


-3.10 


7.13 


0.11 


-2.96 


6.91 


0.11 


-3.14 


-2.58 


0.19 


-4.13 


-2.58 


0.20 


-4.45 


-4.13 


0.26 


as 


a function of Mgb 


-3.37 


0.55 


0.09 


-3.22 


0.54 


0.09 


-3.04 


-0.20 


0.19 


-4.02 


-0.20 


0.20 


-4.29 


-0.32 


0.26 


as 


a function of H7yi 


-2.69 


-0.24 


0.17 


-2.57 


-0.23 


0.16 


-3.33 


0.08 


0.21 


-4.32 


0.08 


0.23 


-4.76 


0.12 


0.31 


as 


a function of C4668 


-2.37 


0.19 


0.15 


-2.25 


0.18 


0.15 


-3.35 


-0.08 


0.16 


-4.34 


-0.08 


0.17 


-4.81 


-0.12 


0.21 



Note. — The fit is M = co + ci x (galaxy property) for models with Z = {0.004, 0.008, 0.02, 0.05} and ages of 5, 8, 12, 17 Gyr. The rms scatter about the fitted 
line is <T, given in magnitudes. 



Table 4 



Linear coefficients for SBF fc-coRRECTiONS as a Function of Model Metallicity (Z) for cz < 15, 000 km s 

























z 




Ic 


F814W 


J 


F160W 


K' 


K 


V-Ic 


V - K 


J -K 


0.004 


3.53 (0.08) 


3.79 (0.30) 


3.24 (0.33) 


-0.05 (0.28) 


0.74 (0.14) 


-2.45 (0.14) 


-2.89 (0.54) 


1.07 (0.02) 


4.32 (0.19) 


2.11 (0.06) 


0.008 


3.62 (0.16) 


6.20 (0.50) 


5.23 (0.52) 


0.33 (0.12) 


1.31 (0.13) 


-2.17 (0.17) 


-2.85 (0.49) 


0.90 (0.03) 


4.59 (0.23) 


2.27 (0.09) 


0.02 


3.40 (0.07) 


6.91 (0.30) 


6.24 (0.20) 


0.25 (0.17) 


1.90 (0.21) 


-0.90 (0.51) 


-2.04 (0.27) 


0.88 (0.06) 


4.88 (0.20) 


2.27 (0.10) 


0.05 


3.65 (0.29) 


6.10 (0.50) 


5.33 (0.55) 


1.02 (0.48) 


2.24 (0.54) 


3.52 (0.45) 


0.61 (0.47) 


1.10 (0.08) 


4.96 (0.29) 


1.97 (0.09) 



Note. — The k(z) term is given by a X z where z is the ledshift and a is the number tabulated above, for different observables (SBF magnitude or galaxy color) as 
a function model metallicity. The convention for the correction is fc(z) = X{z) — X{z = 0) where X is an SBF magnitude or galaxy color. The number in parenthesis 
is the rms scatter for the lineaj: coefficient averaged over models of different ages (see § 3.4). 



Table 5 
Solar Metallicity Models 



Z Gyr 



F%\iW FllOM FWQW J FieOW H 



K' 



K F121M 



M V-Ic V-K J-K H/3 Mg2 Mgb YL^a C4668 



Padova Evolutionary Tracks, Semi-Empirical SEDs {Standard Model) 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 



1 
2 
3 
5 
8 
12 
17 



1.81 
2.36 
2.60 
2.83 
2.61 
2.48 
2.51 



0.46 
0.78 
1.00 
1.25 
1.32 
1.43 
1.52 



-0.77 
-0.35 
-0.05 
0.20 
0.27 
0.40 
0.51 



-2.22 
-1.81 
-1.54 
-1.30 
-1.22 
-1.07 
-0.95 



-2.10 
-1.71 
-1.45 
-1.22 
-1.14 
-1.00 
-0.89 



-3.59 
-3.64 
-3.72 
-3.53 
-3.44 
-3.44 
-3.41 



-3.51 
-3.48 
-3.51 
-3.31 
-3.23 
-3.20 
-3.15 



-3.91 
-3.99 
-4.06 
-3.87 
-3.79 
-3.77 
-3.74 



-5.01 
-5.05 
-5.09 
-4.89 
-4.79 
-4.76 
-4.73 



-5.16 
-5.22 
-5.27 
-5.07 
-4.96 
-4.95 
-4.93 



-5.64 -5.65 -5.69 
-5.78 -5.80 -5.83 
-5.80 -5.82 -5.83 
-5.59 -5.62 -5.61 
-5.47 -5.50 -5.49 
-5.44 -5.48 -5.47 
-5.41 -5.46 -5.45 



-5.64 
-5.82 
-5.88 
-5.68 
-5.56 
-5.58 
-5.58 



-7.09 -7.12 -7.50 
-6.89 -6.92 -7.31 
-6.65 -6.70 -7.05 
-6.20 -6.28 -6.59 
-5.96 -6.05 -6.33 
-5.84 -5.99 -6.26 
-5.80 -5.98 -6.26 



0.82 
1.04 
1.13 
1.17 
1.23 
1.29 
1.33 



2.19 
2.72 
2.94 
2.99 
3.12 
3.25 
3.34 



0.70 
0.84 
0.88 
0.87 
0.90 
0.92 
0.93 



3.86 
2.79 
2.31 
2.04 
1.77 
1.54 
1.37 



0.12 
0.17 
0.20 
0.22 
0.25 
0.27 
0.29 



1.66 
2.50 
3.03 
3.36 
3.72 
4.04 
4.24 



4.42 
-0.30 
-3.16 
-4.52 
-5.70 
-6.71 
-7.44 



2.86 
3.97 
4.73 
5.04 
5.47 
5.85 
6.14 



Padova Evolutionary Tracks, Theoretical SEDs (Difference from Standard Model) 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 



1 
2 
3 
5 
8 
12 
17 



-0.01 
-0.07 
-0.09 
-0.09 
-0.06 
-0.04 
-0.04 



-0.07 
-0.13 
-0.13 
-0.14 
-0.11 
-0.11 
-0.11 



-0.05 
-0.07 
-0.07 
-0.08 
-0.05 
-0.05 
-0.05 



0.07 
0.15 
0.16 
0.17 
0.15 
0.16 
0.17 



0.06 
0.14 
0.16 
0.16 
0.15 
0.16 
0.17 



0.12 
0.09 
0.09 
0.09 
0.07 
0.07 
0.07 



0.08 
0.06 
0.05 
0.05 
0.03 
0.03 
0.02 



0.01 
-0.05 
-0.08 
-0.08 
-0.09 
-0.11 
-0.13 



-0.06 
-0.07 
-0.08 
-0.09 
-0.08 
-0.08 
-0.08 



-0.06 
-0.07 
-0.07 
-0.07 
-0.06 
-0.06 
-0.06 



-0.01 
0.00 
0.00 
0.00 
0.00 
0.00 

-0.01 



-0.01 
0.00 
0.00 
0.00 
0.00 
-0.01 
-0.01 



0.00 
0.01 
0.01 
0.00 
0.00 
-0.01 
-0.01 



0.00 
0.01 
0.01 
0.00 
0.00 
-0.01 
-0.02 



-0.01 
0.01 
0.03 
0.04 
0.05 
0.06 
0.06 



-0.01 
0.01 
0.03 
0.04 
0.05 
0.05 
0.05 



-0.01 0.01 
0.00 -0.03 
0.01 -0.04 
0.02 -0.05 
0.02 -0.06 
0.02 -0.07 
0.02 -0.07 



0.02 
0.01 
0.00 
0.00 
0.00 
0.01 
0.01 



0.01 
0.01 
0.00 
0.00 
0.00 
0.00 
-0.01 



-0.05 
-0.02 
0.00 
0.01 
0.01 
0.01 
0.01 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.05 
0.06 
0.07 
0.06 
0.05 
0.05 
0.05 



-0.15 
-0.14 
-0.10 
-0.08 
-0.06 
-0.04 
-0.01 



0.12 
0.14 
0.15 
0.14 
0.13 
0.12 
0.12 



Padova Evolutionary Tracks, Empirical SEDs (Difference from Standard Model) 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 



1 

2 
3 
5 
8 
12 
17 



-0.04 
-0.11 
-0.13 
-0.12 
-0.06 
-0.14 
-0.18 



0.00 
0.01 
0.01 
0.02 
0.02 
-0.01 
-0.01 



0.03 
0.06 
0.07 
0.07 
0.06 
0.05 
0.05 



0.03 
0.04 
0.07 
0.07 
0.07 
0.08 
0.09 



0.04 
0.06 
0.08 
0.09 
0.10 
0.11 
0.12 



0.06 
0.10 
0.12 
0.13 
0.12 
0.16 
0.18 



0.02 
0.04 
0.05 
0.06 
0.05 
0.06 
0.07 



-0.06 
-0.07 
-0.09 
-0.09 
-0.09 
-0.11 
-0.12 



-0.10 

-0.14 
-0.17 
-0.16 
-0.15 
-0.16 
-0.16 



-0.10 

-0.14 
-0.17 
-0.15 
-0.15 

-0.15 -0.01 
-0.15 -0.01 



0.03 -0.04 -0.04 
0.01 -0.01 -0.01 
0.02 -0.02 
0.00 -0.01 
0.00 -0.01 
0.01 
0.01 



0.01 
0.00 
0.00 
0.00 
0.00 



-0.05 
-0.03 
-0.03 
-0.02 
-0.02 
-0.01 
-0.01 



-0.02 
0.02 
0.05 
0.10 
0.12 
0.13 
0.12 



-0.02 
0.02 
0.09 
0.16 
0.18 
0.24 
0.25 



0.00 

0.08 
0.16 
0.27 
0.33 
0.46 
0.50 



0.04 
0.02 
0.01 
0.01 
0.01 
0.02 
0.02 



0.02 
0.03 
0.02 
0.04 
0.05 
0.04 
0.03 



0.01 
0.01 
0.00 
0.00 
0.01 
0.00 
0.00 



-0.06 
-0.02 
-0.01 
0.02 
0.01 
0.01 
0.00 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.10 
0.11 
0.10 
0.07 
0.06 
0.06 
0.05 



-0.51 
-0.42 
-0.30 
-0.20 
-0.18 
-0.14 
-0.11 



0.19 
0.19 
0.18 
0.10 
0.10 
0.10 
0.08 



Geneva Evolutionary Tracks, Semi-Empirical SEDs (Difference from Standard Model) 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 



1 
2 
3 
5 

8 
12 
17 



-0.10 
0.20 
0.09 
0.09 

-0.04 
0.06 
0.15 



-0.15 
0.22 
0.11 
0.01 

-0.07 
-0.06 
0.05 



-0.08 
0.17 
0.08 

-0.10 

-0.15 
-0.13 
-0.02 



0.16 
0.03 
0.03 
-0.23 

-0.31 
-0.25 
-0.08 



0.14 
0.05 
0.04 
-0.21 

-0.29 
-0.22 
-0.05 



0.47 
0.11 
0.30 
0.09 

-0.10 
0.13 
0.47 



0.41 
0.08 
0.24 
0.02 
-0.16 
0.04 
0.36 



0.42 
0.11 
0.29 
0.07 
-0.12 
0.10 
0.44 



0.41 
0.13 
0.32 
0.12 
-0.07 
0.14 
0.46 



0.41 
0.15 
0.34 
0.15 
-0.04 
0.19 
0.52 



0.38 
0.16 
0.34 
0.19 
-0.02 
0.18 
0.51 



0.38 
0.17 
0.35 
0.21 
-0.01 
0.21 
0.55 



0.35 
0.16 
0.34 
0.21 
0.00 
0.22 
0.55 



0.40 
0.20 
0.40 
0.25 
0.03 
0.28 
0.64 



0.22 
0.12 
0.20 
0.29 
0.04 
0.20 
0.49 



0.22 
0.14 
0.24 
0.35 
0.11 
0.33 
0.64 



0.16 
0.12 
0.20 
0.38 
0.14 
0.34 
0.64 



0.06 
0.00 
0.00 
0.07 

0.09 
0.08 
0.06 



0.07 
0.02 
-0.06 
0.12 

0.16 
0.08 
-0.02 



0.01 
0.01 

-0.02 
0.03 
0.03 
0.00 

-0.03 



-0.01 
0.26 
0.26 
0.14 
0.06 
0.02 
0.00 



0.02 
-0.01 
-0.01 
0.00 
0.01 
0.01 
0.01 



0.01 
-0.18 
-0.23 
-0.01 
0.09 
0.10 
0.11 



0.40 
1.00 
1.33 
0.74 
0.29 
0.07 
0.11 



0.01 
-0.47 
-0.46 
-0.13 
0.09 
0.10 
-0.01 



Geneva Evolutionary Tracks, Theoretical SEDs (Difference from Standard Model) 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 



1 
2 
3 
5 
8 
12 
17 



-0.10 
0.12 
-0.01 
-0.05 
-0.13 
-0.02 
0.08 



-0.16 
0.04 
-0.05 
-0.18 
-0.23 
-0.22 
-0.10 



-0.07 
0.07 
-0.02 
-0.20 
-0.23 
-0.22 
-0.11 



0.19 
0.13 
0.15 
-0.15 
-0.23 
-0.19 
-0.03 



0.17 
0.14 
0.16 
-0.13 
-0.21 
-0.16 
0.01 



0.49 
0.23 
0.41 
0.17 
-0.04 
0.19 
0.53 



0.42 
0.16 
0.31 
0.06 
-0.13 
0.07 
0.39 



0.42 
0.11 
0.28 
0.05 
-0.15 
0.07 
0.41 



0.41 
0.07 
0.24 
0.06 
-0.13 
0.09 
0.41 



0.42 
0.09 
0.28 
0.10 

-0.09 
0.14 



0.40 
0.18 
0.35 
0.22 

0.00 
0.21 



0.39 
0.20 
0.37 
0.24 
0.02 
0.25 



0.36 
0.19 
0.36 
0.25 
0.03 
0.25 



0.47 0.54 0.58 0.58 



0.41 
0.23 
0.43 
0.29 
0.07 
0.32 
0.68 



0.22 
0.14 
0.22 
0.34 
0.11 
0.28 



0.22 
0.15 
0.26 
0.41 
0.17 
0.40 



0.16 
0.13 
0.21 
0.42 
0.18 
0.39 



0.57 0.72 0.70 



0.06 
-0.04 
-0.05 
0.00 
0.01 
0.00 
-0.03 



0.11 
0.02 

-0.07 
0.10 
0.14 
0.06 

-0.03 



0.02 
0.02 

-0.01 
0.03 
0.03 
0.00 

-0.03 



-0.07 
0.24 
0.25 
0.15 
0.08 
0.05 
0.04 



0.02 
0.00 
-0.01 
0.01 
0.01 
0.01 
0.01 



0.05 
-0.06 
-0.12 
0.11 
0.21 
0.23 
0.24 



0.21 
0.86 
1.21 
0.65 
0.23 
0.05 
0.12 



0.13 
-0.26 
-0.25 
0.10 
0.31 
0.32 
0.21 



Geneva Evolutionary Tracks, Empirical SEDs (Difference from Standard Model) 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 



1 
2 

3 
5 
8 
12 
17 



-0.15 
0.14 
-0.02 
-0.05 
-0.16 
-0.13 
-0.04 



-0.14 
0.24 
0.11 
0.01 
-0.08 
-0.07 
0.04 



-0.03 
0.23 
0.11 
-0.06 
-0.12 
-0.09 
0.02 



0.18 
0.06 
0.03 
-0.20 
-0.28 
-0.20 
-0.04 



0.18 

0.08 
0.06 
-0.17 
-0.25 
-0.16 
0.01 



0.46 
0.17 
0.37 
0.14 
-0.04 
0.18 
0.52 



0.39 
0.11 

0.28 
0.05 
-0.13 
0.07 
0.38 



0.38 
0.06 

0.25 
0.02 
-0.17 
0.04 
0.36 



0.39 
0.02 

0.18 
-0.02 
-0.22 
-0.02 

0.31 



0.39 
0.03 

0.21 
0.01 
-0.19 
0.02 
0.36 



0.41 

0.16 
0.34 
0.19 
-0.02 
0.18 
0.51 



0.40 
0.17 

0.36 
0.21 
0.00 
0.20 
0.53 



0.37 
0.16 

0.35 
0.22 
0.01 
0.21 
0.54 



0.41 
0.18 
0.40 
0.24 
0.03 
0.26 
0.61 



0.22 
0.14 
0.25 
0.42 
0.21 
0.37 
0.64 



0.21 
0.15 
0.30 
0.50 
0.30 
0.52 
0.81 



0.21 
0.19 
0.31 
0.62 
0.46 
0.70 
0.97 



0.07 
0.02 
0.01 
0.08 
0.10 
0.07 
0.06 



0.15 
0.07 
-0.02 
0.17 
0.21 
0.11 
0.01 



0.05 
0.03 

-0.01 
0.03 
0.04 
0.01 

-0.02 



-0.08 
0.27 
0.26 
0.15 
0.07 
0.04 
0.02 



0.02 

-0.01 
-0.01 
0.00 
0.01 
0.01 
0.01 



0.13 

-0.04 
-0.08 
0.12 
0.21 
0.19 
0.19 



-0.19 
0.73 
1.03 
0.49 
0.07 

-0.02 
0.05 



0.24 
-0.31 
-0.26 
0.07 
0.27 
0.22 
0.09 



Note. — The top set of results gives the model output lor our standard solar-metallicity model (Padova evolutionary tracks with semi-empirical SEDs). For all other models, the differences 
from our standard models are tabulated, e.g., {Geneva model. Pickles SED} — {Padova model, semi-empirical SEDs}. 



